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Abstract 


This  paper  Is  the  first  In  a  series  of  three  In  which  we  dlscuss^ome 
theoretical  and  practical  aspect  of  a  feedback  finite  element  method  for 
solving  systems  of  linear  second  order  elliptic  partial  differential  equations 
(with  particular  interest  in  classical  linear  elasticity) .  ^^n^this  first  part 
^-we-  Introduce ipame  nonstandard  finite  element  spaces,  though  based  on  the  usual 
square  bilinear  elements,  permit  local  mesh  refinement.  The  algebraic 


structure  of  these  spaces  and  their  approximation  properties  are  analysed.  An 
^ ""equivalent  estimator^for  the  H*,  finite  element  error  is  developed.  ^  In  the 
second  paper  we  shall  discuss  the  asymptotic  properties  of  this  estimator.  In 
the  third  paper  we  shall  also  report  on  some  computational  experience  with  the 
FEARS  program  which  uses  this  estimator  as  part  of  a  feedback  loop  to  control 
mesh  refinement  and  some  of  its  programing  features. 


§0.  Introduction 

The  practical  success  or  failure  of  many  finite  element  computations 
often  depends  critically  on  the  user's  choices  of  finite  element  mesh  and 
element  type.  As  a  simple  illustration  of  this,  consider  the  boundary  value 
problems  that  arise  in  classical  plane  linear  elasticity.  For  such  problems 
it  is  well  known  that  in  the  neighbourhood  of  certain  crltial  boundary  points 
(e.g.  angular  boundary  points,  or  points  where  the  boundary  conditions  change 
between  specified  tractions  and  specified  displacements),  the  stresses  exhibit 
some  form  of  singular  behaviour.  Unless  such  critical  points  are  handled 
carefully,  the  resulting  finite  element  solution  may  have  disappointing  ac¬ 
curacy  . 

One  way  in  which  such  critical  points  may  be  treated  is  to  employ  an 
appropriate  mesh  refinement  strategy  in  the  neighbourhood  of  the  point. 

Broadly  speaking,  two  kinds  of  refinement  strategy  can  be  identified.  On  the 
one  hand,  there  are  a  priori  refinement  techniques  which  grade  the  mesh  in  a 
manner  governed  either  by  earlier  experience  with  similar  kinds  of  problems  or 
by  the  results  of  some  a  priori  analysis  of  the  nature  of  the  singularity. 

For  many  problems  in  linear  elasticity  an  asymptotic  representation  of  the 
solution  in  the  vicinity  of  the  critical  point  is  available.  Using  such 
representations  it  is  often  possible  to  derive  a  sequence  of  graded  meshes 
which  can  be  shown  to  converge  in  the  energy  norm  at  an  optimal  rate  with 
respect  to  the  number  of  degrees-of-freedom  of  the  resulting  discrete  system. 
From  a  practical  point  of  view  however,  a  weakness  of  such  a  priori  methods  is 
that  the  analysis  or  experience  they  are  based  upon  is  asymptotic  in  nature 
and  is  seldom  discriminating  enough  to  tell  whether  a  singularity,  though 
present  in  theory,  is  going  to  cause  significant  problems  at  the  level  of 
accuracy  that  one  is  working.  As  an  example  consider  the  stress  singularity 
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Kr~  ^  g(6)  typically  associated  with  cracks  in  plane  elasticity.  The 
potential  of  this  singularity  to  affect  the  accuracy  of  a  finite  element 
approximation  will  depend  on  K.  Usually  an  a  priori  analysis  can  give  little 
insight  into  the  value  of  K.  If  K  is  small  enough  (compared  to  the  overall 
level  of  accuracy  desired),  then  no  harm  will  be  done  by  employing,  say,  a 
uniform  mesh  near  the  crack  tip.  Indeed,  were  a  refined  mesh  used,  the  extra 
degrees  of  freedom  would  not  lead  to  any  significant  Improvement  in  accuracy 
over  the  uniform  case.  On  the  other  hand,  any  (fixed)  mesh  refined  near  the 
crack  tip  will  give  a  far  from  optimal  mesh  as  K  ®.  (Optimal  here 
indicates  a  sense  of  minimum  error  for  a  given  number  of  degrees  of 
freedom).  This  kind  of  phenomenon  becomes  more  pronounced  when,  as  occurs  in 
most  practical  problems,  there  are  a  number  of  critical  boundary  points  in  the 
region  of  Interest,  and  a  decision  must  be  made  on  how  the  corresponding 
refinements  should  be  "weighted." 

The  other  kind  of  refinement  strategy  referred  to  earlier  is  a  posteriori 
in  character.  In  this  approach  an  initial  finite  element  solution  is 
calculated  using  some  mesh.  This  solution  is  then  examined  in  some  fashion, 
and  based  upon  this  examination  a  refinement  of  the  initial  mesh  is  decided 
upon.  Using  this  new  mesh,  a  new  finite  element  solution  is  computed.  This 
process  can  obviously  be  iterated  until  some  stopping  criterion  is 
satisfied.  This  kind  of  feedback  technique  could  conceivably  avoid  the 
problems  mentioned  above,  since  the  computational  significance  of  the 
singularities  present  may  be  able  to  be  ascertained  during  the  feedback 
process.  Everything,  of  course,  depends  on  the  a  posteriori  examination 
carried  out  after  each  step  and  the  refinement  decision  arising  from  it. 

So  far  we  have  only  mentioned  mesh  refinement  necessitated  by  some  form 
of  singular  behaviour  of  the  solution.  However  in  many  other  situations 
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proper  mesh  refinement  can  also  be  crucial.  As  a  typical  example  consider  the 
use  of  curved  elements  for  solving  problems  In  elasticity.  In  general  such 
elements  cannot  exactly  represent  rigid  body  motion.  Because  of  this,  proper 
mesh  refinement  may  again  be  essential  for  satisfactory  results,  even  though 
the  solution  Is  very  smooth.  It  would  seem  very  difficult  to  predict  the 
pattern  of  such  refinement  a  priori. 

In  this  series  of  papers  we  shall  describe  and  provide  an  analysis  of  one 
such  feedback  approach.  The  approach  we  shall  deal  with  has  been  Implemented 
in  a  practical  form  in  the  FEARS  program  [see  [1]}.  This  algorithm  is  based 
upon  an  a  posteriori  examination  which  involves  the  calculation  of  an  a 
posteriori  error  estimator  for  the  energy  (or  similar)  norm  of  the  error. 

This  estimator  is  composed  of  elementwise  error  "indicators",  and  the 
refinement  decision  at  each  stage  is  made  on  the  basis  of  the  distribution  of 
these  indicators.  The  theoretical  analysis  of  this  method  is  far  from 
complete.  At  the  moment  there  are  many  cojectures,  etc.,  which  though 
convincingly  demonstrated  by  many  numerical  experiments  can  either  not  be 
proved  rigorously,  or  not  be  proved  in  the  generality  that  practical 
experience  would  indicate  they  hold. 

Our  theoretical  analysis  will  concentrate  upon  the  error  indicators  and 
estimators,  in  particular,  upon  clarifying  in  what  sense  they  "estimate"  the 
energy  norm  of  the  error  of  the  finite  element  solution.  We  shall  show  under 
quite  weak  assumptions  that  the  estimator,  E  say,  is  an  equivalent  estimator 
for  the  energy  norm  lei  of  the  error  e,  that  is 

c“lE  <  lei  <  CE  (1) 

for  some  constant  C  >  0  uniformly  over  large  class  of  meshes;  and  moreover 
under  some  further  restrictions  that  E  is  also  an  asymtotically  exact  esti¬ 
mator  for  lei,  that  is 
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lei  •  E(l+o(lel))  (2) 

as  lei  ♦  0. 

The  properties  (1)  and  especially  (2)  suggest  that  the  estimator  can  be 
reliably  used  for  stopping  the  refinement  process  once  some  desired  accuracy 
Is  achieved.  Experience  with  FEARS  confirms  this. 

To  be  considered  worthwhile  the  process  of  succesive  construction  of 
meshes  should  lead  to  a  sequence  of  meshes  whose  rate  of  convergence  Is 
comparable  to  that  of  the  theoretically  optimal  mesh  grading  (at  least  when 
dealing  with  practical  problems).  They  should  also  have  other  "optimal" 
properties.  When  a  feedback  process  has  such  optimal  properties,  it  is  called 
an  adaptive  process.  (For  more  about  this  see  [2]  [3]  [4].)  Experience 
indicates  that  the  FEARS  program  implements  an  adaptive  process. 

In  §1.1  of  this  paper  we  shall  describe  the  kinds  of  boundary  value 
problems  that  we  wish  to  consider  (essentially  those  related  to  linear  elas¬ 
ticity).  §1.2  contains  a  description  of  the  finite  element  discretization  to 
be  employed.  In  §1.3-1. 5  we  give  some  properties  of  the  corresponding  finite 
element  spaces  for  use  later  in  the  paper.  In  particular,  we  introduce  the 
important  concept  of  a  K-mesh  in  §1.4.  Although  we  describe  a  rather  general 
set  up  in  §1.1  and  §1.2,  the  subsequent  analysis  is  carried  out  in  a  more 
restricted  setting.  This  has  been  done  for  the  sake  of  clarity  and  simplicity 
of  notation.  The  analysis  of  the  general  case  can  be  done  analogously.  In  §2 
we  derive  an  equivalent  estimator  for  the  energy  norm  of  the  finite  element 
error.  An  important  step  in  this  is  the  basic  error  estimate  of  §2.1.  Using 
this  result  the  error  is  able  to  be  localized  to  a  small  number  of  elements. 
Some  technical  lemmas  are  proved  in  §2.2  and  §2.3,  while  §2.4  contains  the 
main  results  of  this  section. 

In  the  second  paper  of  this  series  we  shall  deal  with  the  asymptotic 
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exactness  of  the  estimator,  as  well  as  discussing  the  overall  performance  of 
the  algorithm. 

The  third  part  will  deal  more  specifically  with  design  of  the  FEARS 
program  and  analyse  its  performance  in  the  light  of  the  developed  theory. 


Formulation  of  the  problem  and  its  finite  element  solution 


1.1.  We  shall  consider  the  boundary  value  problem 


Li(u)  = 


2  2 

'j.Ji-i  +  £  ci3“j  '  fi  ln  0 


u^  ■  0  on  a  Q, 


2  2 


Jj  '’ik<u)“k  5  Jj  aiju0iujK  ■  ci  °n  al°  (1  -  >. 

on  a  bounded  domain  Q  c  Jr  whose  boundary  3Q  is  made  up  of  two  disjoint 
parts  J4  0  and  9*Q;  n  ■  (nj,^)  is  the  outward  pointing  unit  normal 

on  9Q.  This  problem  can  be  cast  in  a  variational  form.  If  d^Q  is 
sufficiently  smooth  (piecewise  smooth,  say)  we  may  define  the  trace  of  H1 (Q) 
functions  on  Let  us  write 

H  -  {(Vj^)!  v±  €  H^Q)} 

Hq  ■  {(vltv2):  v±  €  H*(Q),  -  0  on  3°Q}. 

and  Hq  are  Hilbert  spaces  with  respect  to  the  norm 
2  2  1/ 

lvlL  q  ■  (  £  *vi’l  q)  ^  »  with  *vi’l  q  being  the  usual  (scalar)  H1(Q) 

Sobolev  norm.  The  boundary  value  problem  (1.1)  can  now  be  posed  as:  Find 
u  €  H  o  such  that 

2  2 

b(u,v)  -  /  l  f.v,  d x  +  /  l  t±v.  ds,  Vv  €  HQ  (1. 

01-1  5Xo  1-1 

where  the  bilinear  form  b:  Hq  x  Hq  -*■  R  is  defined  by 


*  L 

b(w,z)  -  /(  l  aiiklDiw1Dkzi  +  X  ci1w1zi)dx* 

Q  i,j,k,i-l  *  3  K  1  i,j-l  1J  J  1 


We  shall  assume  that  Q  can  be  represented  as  a  collection  of 

transformed  unit  squares.  To  this  end  we  make  the  following  assumption  on  Q: 

There  Is  a  finite  number  of  subdomains  Q.  £  Q(d  *1 . N)  such  that 

_____  £ 

(1)  QJfl  -  0  (d,g  -  1 . .  d  i  g) 

d  g 

N 

(11)  Q  -  U  Q 
d-1  d 

(111)  There  Is  an  invertible  transformation  4^:  [0,1]^  which, 

together  with  is  Inverse  41,^,  is  sufficiently  smooth  (bounded 

a 

derivatives  of  all  orders,  say) 

(lv)  Qd  -  ^((O.l)2). 

— j  2 

The  image  under  4»j  of  the  closed  edges  and  corners  of  [0,1]  will  be 
referred  to  as  the  edges  and  corners  respectively  of  Q^. 

(v)  Qd  and  Qe(d  j4  e)  can  only  have  a  single  edge  or  a  single  corner 

In  common  (or  else  Q.  and  Q  are  disjoint) 

a  e 

(vl)  d^Q  is  the  union  of  a  number  of  (complete)  edges  of  subdomains. 

We  shall  write  9?(0,1)2  -  «|»,O0Q  fl  QQ  )  and  aj(0,l)2 
a  a  a  a 

-  4»d(aio  n  e>Qd). 

An  edge  common  to  Q .  and  Q  will  be  called  an  interface  between  Q .  and  Q 

d  e  ■  ■  ■  d  e 

(vll)  If  r  is  an  interface  between  Q.  and  Q  ,  then  4>,  and  4; 

a  e  a  e 

must  "agree”  on  T  in  the  following  sense:  if  s  is  an  arclength 
measured  along  T,  and  sd  and  se  are  arclengths  measured  along 
4>d(r)  and  4<e(r)  respectively,  then  regarding  sd  and  se  as 
functions  of  s,  either  sd(s)  ■  se(s)  or  sd(s)  ■  1  -  se(s). 

Figure  1  shows  a  possible  partitioning  of  a  circle  into  subdomains.  For 
a  suitable  choice  for  the  mappings  4»,  in  the  case,  see  [5].  If  f  is  an 


Interface  between  Q.  and  Q  ,  then  we  can  naturally  Identify  Q.(r)  and  Q  (F) 

d  e  d  e 

as  illustrated  in  Fig.  1.  By  virtue  of  (vii)  above  this  identification 
takes  a  particularly  simple  form  when  expressed  in  terms  of  local  arclength. 


Figure  1.  Partitioning  of  the  circle  into  subdomains. 

This  representation  of  Q  induces  a  natural  correspondence  between  H 
and  the  set  M  of  N  tuples  of  the  form 

V  -  (V^1),...,V^d),...,V(N)) 

where 

(i)  V(d)  -  (v[d),V^d)),  V^d)  €  H^CO.l)2)  (i  -  1,2,  d  -  1 . N) 

(ii)  if  T  is  an  interface  between  Q,  and  Q  (d,e  -  1,...,N;  di*e) 

d  e 


Vj  o  *d  -  Vj  °  ♦.  0n  r 


(J  “  1,2) 


This  correspondence  is  defined  by: 

v  “  (vl »v2>  €  H  -*•  V  -  (....(Vj  o  v2  o  t^1),...)  6  M, 


W  -  (W(1\...,W(N))  f  I  +  w  -  (wj.w 2)  €  tf 


where 


wi  "  wi  ^  °  on  °d’  (1  “  1,2»  d  “  1»***»N> 

(w^ €  H*(Q)  by  virtue  of  the  Interface  continuity  requirement  (11)  above). 

Let  Mq  £  M  be  the  set  of  tuples  satisfying  the  additional  condition 
(ill)  If  T  is  an  edge  of  Q.  and  T  £  B^Q  then 


v‘d>  o^-O  on  r. 


(j  -  1,2). 


Clearly  Mq  •«-*>  Hq  under  the  above  natural  correspondence.  M  and  Mq  are 
Hilbert  spaces  with  respect  the  norm 


N  2 


'l  fQ  q  ]  “  (I  I  »^d)l2  2f2- 

d-1  i-i  1  1,(0,1; 


The  problem  (1.2)  can  now  be  reformulated  as:  Find  U  €  Mq  such  that 
B(U,V)  -  f  (  /  l  F<d)V<d)  dx  +  /  l  T<dVd)  ds),  VV  (H 

d-L  (o.i)21-1  ako.D21-1 

d  (1.3a) 

where  the  bilinear  form  B:  Mq  *  Mq  R  is  defined  by 


B(W,Z)  -  l  I  (  l  A^klDXW1d)Vld)+  2  C^Va'zJ°')dx.  (1. 
d-1  ,n  ,,2  i.j.k.l-l  1  J  *  i,j«l  J 


,(d)„(d)_(d)i 


In  (1.3)  we  have  used  the  notation 


Lijia 

|E(d)| 

r(d)  . 

1 

■  (C1J  °  *dl)- 

Cij 

|E(d)| 

r(d)  . 

1 

(ft  O  4,"1), 

Fi 

|E(d)| 

(0,l)2; 

while 

on  &J(0,1) 

II 

/—V 

T* 

1 

E(d) 

(ti  °  ^dl>* 

We  have  set  -  Di^d^j*  lE^  I  *  determinant  (e£^:  i,j  ■  1,2),  and 

f(d)  »  (  £  ((J»d)i)^)/^  with  denoting  differentiation  with  respect 

1  2 

to  arclength  along  a  (0,1)  . 

U 

With  regard  to  the  coefficients  and  Input  data  of  (1.2)  we  shall  suppose 
that  the  cij  and  are  sufficiently  smooth  on  each  separately 

(say,  bounded  derivatives  of  all  orders),  and  that  the  t^  are  sufficiently 
smooth  on  each  subdomain  edge  contained  in  d^Q. 

We  assume  the  symmetries 

^4  41,1  m  a44  01*  (1|  Jl^tl  m  1»2), 


so  ensuring  that  b(»,«),  is  a  symmetric  bilinear  form.  Additionally,  we 
shall  suppose  that  the  bilinear  form  b(*,*)  is  coercive  over  Hq,  that  is, 
there  exists  a  >  0  such  that 


(1.4a 


|b(w,w) |  >  alwl1  Q  iw  (  H( 


Further  we  will  assume  that  for  some  C  >  0 


2 

VjkkMj  >  ««1  +  «2>>  k  ■  ».»!  5  «  *2. 


(1.4b 


By  virtue  of  our  assumptions  on  the  mappings  <1^,  the  above  properties 
transfer  naturally  to  the  transformed  system.  Let  us  explicitly  note  the 
symmetry  condition  B(0,V)  -  B(V,U)  VV,W  €  Mq,  the  coercivity  condition 


|B(W,W)  |  >  «IWI,  VW€M0, 


t }jml 


k  -  1,2;  ?  k 


Mote  that  In  the  case  when  the  a.  ...  are  discontinuous  across  an 

ljKX 

Interface  T  between  two  subdomains,  then  the  classical  formulation  (1.1) 
needs  to  be  supplemented  by  an  Interface  condition  expressing  the  continuity 
of  "tractions"  across  T.  This  condition  Is  of  course  implicit  In  the  varia¬ 
tional  formulations  (1.2)  and  (1.3)  of  the  problem. 

The  finite  element  approximation  that  we  shall  discuss  Is  based  upon  the 
formulation  (1.3)  of  the  problem.  If  c  MQ,  the  corresponding  Galerkln 
approximation  U  €  Mq  to  U  is  defined  by 


N  2  2 

B(0,V)  -  I  (  /  |  F(d)V(d)dx  +  /  l  Tjd)V.(d)ds)  W  €  Hq 

J _ 1  e  j  11  m  t  t  1  1  w 


s 


The  coercivity  of  B(» ,  • )  ensures  Che  existence  and  uniqueness  of  tJ.  We 
have  of  course  Che  projection  property 


lW,i,(fllt...,od)  *  c,uli,(Q1,...,Qdr 


(1.5 


The  finite  element  error  U-tf  satisfies  the  usual  orthogonality  relation 


B(U-tf,V) 


VV€EL 


(1.6a 


which  leads  to  the  standard  kind  of  best  approximation  estimate 


'^'1.(0. . a.)  ‘  Clnf  ,u"vli,(a. . a.)- 

1  d  vefl0  1  d  . 


The  finite  element  approximation  U  corresponds  (under  the  natural 
correspondence  between  H  and  M)  to  an  approximation  u  to  the  solution 
u  of  (1.1),  (1.2).  In  terms  of  the  energy  norms  of  the  respective  errors  we 
have 

(b(u-u,u-u))/2  -  (B(U-U,U-U))1/z  .  (1.7 


1.2.  Suppose  that  a  subdivision  of  Q  into  subdomains  Q^,...,QN  with 
corresponding  mappings  <J^,...,4»N  satisfying  the  conditions  of  §1.1  has  been 
decided  upon.  We  shall  now  define  what  we  mean  by  a  mesh  p  on  Q. 

A  jie8h_2_  is  an  N-tuple  (Pj,...,Pjj)  where  each  Pd  (d  ■  1,...,N)  is  a 
partition  of  [0,1]4  into  closed  squares  A  (with  edges  parallel  to  the 
coordinate  axes).  Each  Pd  is  called  a  submesh  on  Qd  and  either 

(i)  vd  -pg-  {[0,1]2} 

or 

(ii)  P<j  ■  P^^is  constructed  from  an  existing  submesh 

P£*“^t  by  replacing  any  A  €  P^*  ^  by  the  four  congruent 


squares  resulting  froa  the  simultaneous  bisections  of  A  In  the 
two  coordinate  directions  (see  Fig.  2). 


-  (tO.lJ2}  Pj15  -  -  <A2,...,A8} 

Figure  2.  The  construction  of  a  submesh  on  the  subdomain  • 

Each  dosed  square  A  €  P4  Is  called  an  element  of  P4.  He  shall  use 
I A I  to  denote  the  length  of  a  side  of  A.  Clearly  | A I  “  2” 8  with  s  >  0  an 
Integer.  Further  we  denote  h(0,)  -  max  t  A I • 

A  point  P  Is  called  a  node  of  P d  If  either 

(1)  P  Is  a  vertex  of  an  element  of  P4 
or 

(11)  there  Is  an  Interface  between  and  Qe,  and  <j>e  o  <J»^(P)  is  a 

vertex  of  an  element  of  Pe. 

Nodes  P  of  P4  are  classified  as  P-proper  if  either 

(1)  P  €  (0,1)2  and  whenever  P  €  A'  for  A'  €  P<j,  then  P  Is  a 
vertex  of  A' 
or 

(11)  ^(P)  €  dQ 

or 

(111)  4>.*(P)  lies  on  an  Interface  between  Q.  and  Q  ,  P  Is  the  vertex 

a  a  e 

of  an  element  of  P4,  and  (4>e  o  4^  )(P)  Is  the  vertex  of  an 


element  of  p_ 


element  of  Ve . 

If  a  node  is  not  P-proper,  it  is  said  to  be  V -improper.  The  cases  (i),  (ii), 
(ill)  above  are  clearly  mutually  exclusive,  and  we  shall  further  classify  the 
V  -proper  nodes  as  Interior,  boundary  or  interface  nodes  depending  upon 
whether  (1),  (ii)  or  (Hi)  applies.  Fig.  3  shows  an  example  of  mesh  with 
proper  and  Improper  nodes  indicated. 

Each  of  the  straight  line  segments  y 


Figure  3.  A  mesh  on  0  (*-*■  denotes  the  natural  identification  of 

points  across  interfaces). 

forming  the  boundary  of  an  element  of  Pj  will  be  called  an  edge.  For 
definiteness  an  edge  will  be  assumed  to  be  closed,  that  is,  it  Includes  its 
endpoints.  |y|  will  denote  the  length  of  an  edge.  An  edge  is  called  a 
primitive  edge  if  it  contains  no  nodes  other  than  its  endpoints. 

We  associate  with  a  mesh  V  the  finite  element  subspace  W(P) 


consisting  of  all  tuples  of  the  form 


(i)  W(d)  -  (w{d),w£d)),  tf$d)  €  H1<<0,1)2)  (j-1,2;  d-l,...,N) 

(ii)  is  bilinear  on  each  element  A  €  (j-1,2;  d*l,...,N) 

(iii)  if  T  is  an  interface  between  Q,  and  Q 

a  e 

Wjd)  o  <pd  -  W<e)  o  ^  on  r,  (j-1,2;  d,e-l . N). 

Clearly  U(V)  $.  M.  Define  Mq(0)  -  M(V)  (1  Mq  . 

The  finite  element  subspaces  we  have  introduced  are  non-standard  since 
improper  nodes  are  permitted.  (Note  however  that  the  spaces  remain 
conforming).  These  spaces  permit  local  mesh  refinement,  yet  maintain  many  of 
the  desirable  programming  characteristics  of  the  more  usual  square  or 
rectangular  elements. 

Having  defined  a  general  framework,  we  shall  in  the  analysis  for  the 
remainder  of  this  paper  only  consider  the  particular  case  of  one  subdomain 
(N-l)  and  <J>  the  identity  mapping  (so  Q  -  (0,1)^).  This  restriction  is  for 
notational  simplicity  only.  Our  results  extend  quite  naturally  to  the  case  N 
>  1  and  4»  general,  though  with  a  considerable  growth  in  notation.  In  the 
light  of  this  simplification  we  shall  from  now  on  suppress  the  subdomain  index 
and  just  write  Q,  M,  6  (0,1)  ,  etc.  There  is  also  no  need  now  to 
distinguish  between  the  original  problem  (1.2)  and  the  transformed  problem 
(1.3).  For  definiteness  we  shall  from  now  on  use  the  notation  of  the  original 
problem. 

To  further  contain  notation,  where  no  confusion  is  possible  we  shall  not 


notatlonally  distinguish  between  vector  valued  functions  and  their  compo¬ 
nents.  In  such  instances  all  operations,  relations,  etc.,  are  to  be 
understood  in  a  componentwise  sense. 


1.3.  On  a  natter  of  terminology  we  shell  say  that  a  line  [Pj.Pj]  is  a 
bleary  segment  of  another  line  If 


Pl  -  Qi  +  ^  (02-91) 

Ui.1 

p2  -  o,  +  (o2-0!). 


s  *  0,1, ...J  1c  •  0.1.. ... 2^- 


Lewma  3.1. 

(a)  If  A,  A*  are  distinct  elements  of  P  with  |AI  <  |A*|,  then  one 
and  only  one  of  the  following  holds: 


(1)  A  n  A*  -  0 

(11)  A  and  A  share  only  one  point.  This  point  being  a  proper 
node  (and  hence  also  a  vertex  of  both  A  and  A*)  . 


(ill)  A  fl  A*  Is  an  edge  of  A.  This  edge  is  a  binary  segment  of 
an  edge  of  A*. 


(b)  If  Q  Is  an  Improper  node,  then  Q  lies  In  the  Interior  (i.e..  Is 
not  an  endpoint)  of  an  edge  of  a  unique  element  A*.  Furthermore,  If  Q  is 
the  vertex  of  A  €  P  then  | A I  <  I A* I . 


Proof .  (a)  These  results  follow  readily  by  an  induction  based  on  the 

refinement  process  used  to  construct  P. 

(b)  Since  Q  is  Improper,  there  is  an  element  A*  with  Q  €  A  *  but  Q 
is  not  a  vertex  of  *.  It  is  readily  seen  that  *  is  unique.  Case  (iii) 
of  (a)  applies.  Assuming  that  I  I  1  *  1  leads  to  a  contradiction  since 
Q  is  not  a  vertex  of  *.  □ 
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Lemma  3.2.  Suppose  V  is  a  mesh  and  Chat  z  L2((0,l)z)  Is  bilinear  on 
each  A  €  V.  If  z  has  a  well  defined  limiting  value  at  each  Interior  node 
(i.e.,  z  Is  continuous  In  these  points)  then  z  £  $(P). 


Proof.  It  suffices  to  show  that  z  €  C°((0,1)^).  Let  x  6  (0,1)^.  If  x 
lies  In  the  Interior  of  an  element  then  z  Is  obviously  continuous  at  x. 

The  other  possibility  Is  that  x  is  a  common  boundary  point  of  (at  least)  two 
elements,  A  and  A*  say.  There  are  two  cases  to  consider  here:  (i)  x  is 
a  vertex  of  A  or  A*.  In  this  case  the  hypothesis  guarantees  continuity. 
(11)  x  is  not  a  vertex  of  A  or  A *.  Case  (ill)  of  Lemma  3.1(a)  must  then 
apply.  Quite  generally  we  may  suppose  that  |A|  <  |A*I,  so  A*  A  Is  an 
edge,  [Pj , P2 1  say,  of  A  and  is  contained  in  an  edge  of  A*.  The  limit 
of  z  on  [Pj,P2l  from  within  A  is  a  linear  function,  as  is  the  limit  from 
within  A*.  But,  by  our  hypothesis,  at  P^  and  P2  these  limits  must 
agree.  Thus  they  must  agree  througout  and  80  in  particular  at 

x.  The  continuity  of  z  at  x  follows.  □ 


Theorem  3.3. 

(a)  For  any  proper  node  P  there  is  a  <|>p€$(P)  satisfying 


*P(Q) 


1  if  Q  -  P 

0  If  Q  is  a  proper  node,  Q  ^  P. 


(3.1) 


In  addition  <J>p  >  0  everywhere. 

(b)  If  <p  €  S((P)  and  <p(P)  “  0  for  all  proper  nodes  P,  then  <p  ■  0 

(In  particular,  the  are  unique). 

(c)  {4>p!  P  a  proper  node}  is  a  basis  for  W(P).  In  fact,  for  any 

q>  €  M(D),  <p  -  I  <p(PHp  where  V  denotes  summation  over  all  proper  nodes 
P  P 


P  of  V 


Proof,  (a)  We  shall  construct  a  function  z  6  L2(0,l)2  which  is  biliear  on 
each  A  €  P.  Arrange  the  elements  of  D  in  order  A^,...,Am  where 
I A  j  I  >  I  (J  ■  l,...,m-l).  Suppose  that  z  >  0  has  been  defined  on  each 

Aj  for  j  <  n.  To  define  z  on  An  it  suffices  to  specify  the  limiting 
values  of  z  from  withing  Afl  at  each  vertex  Q  of  An: 

(i)  If  Q  is  a  proper  node,  then  set 


z(Q) 


1  if  Q  -  P 
0  if  Q  +  P 


(ii)  If  Q  is  improper,  then  by  Lemma  3.1(b),  Q  lies  in  the 

interior  of  an  edge,  say,  of  some  unique  element  A* 

with  | A* I  >  IAnl.  Thus  A*  ■  Aj  for  some  j  <  n  and  we  may 
define  z(Q)  ■  z*<Q)  >  0  when  z*(Q)  is  the  limiting  value 
of  z  at  Q  from  within  A*. 

Thus  we  are  able  to  define  z  >  0  on  A„.  It  is  clear  from  this  method 

n 

of  construction  that  z  has  a  well  defined  limiting  value  at  each  node  of 
V,  and  so  by  Lemma  3.2,  z  €  W(D).  Obviously  if  we  set  <i>p  *  z,  then 
(pp  >  0  and  (3.1)  is  satisfied. 

(b)  Again  arrange  the  elements  of  V  in  order  of  nonincreasing  size. 
Suppose  that  cp  ■  0  on  all  A;  for  j  <  n.  Consider  any  vertex  Q  of  An. 

If  Q  is  proper  then  <p(Q)  ■  0.  On  the  other  hand,  if  Q  improper,  then 

again  by  Lemma  3.1(b),  Q  lies  in  the  interior  of  an  edge,  [ Pj , P2 1  say*  of 
some  unique  element  A*  with  |A*|  >  IAnl«  Thus  A*  ■  Aj  with  j  <  n. 

But  ip  ■  0,  on  Aj,  and  so  <p(Q)  ■  0.  Thus  'P  ■  0  at  all  vertices  of 

A  ,  and  therefore  <p  ■  0  on  A. 


(c)  This  follows  readily  from  (a)  and  (b) 


2 


Corollary  3.4. 

(a)  {4»p:  P  i  a°(0,l)2}  is  a  basis  for  MQ(P). 

(b)  l  <k  -  1,  14-pl  <  1. 

P 

(c)  On  any  A  €  V,  'Dj^p*  <  (1  “  1*2). 

Proof.  This  follows  readily  from  the  Theorem  3. 

Theorem  3.5. 

2 

(a)  Supp(4>p)  *  closure  {x  €  [0,1]  :  <J>p(x)  f  0}  is  the  union  of  a 
(whole)  number  of  (closed)  elements. 

(b)  Supp(<J»p)  is  "connected"  in  the  sense  that  if  A',  A"  S  supp  <J>p 
then  there  exists  a  sequence  of  elements  A'  ■  Aq,...,As  "  A”  such  that: 

(i)  Aj  £  supp  <J*p  (j  -  0,...,s) 

(ii)  The  pair  A^,  A^  share  an  edge  of  the  smaller  of  the  pair 
(j  ■  0, . . . ,s-l ) . 

Proof,  (a)  This  follows  from  the  simple  observation  that  if  a  bilinear 

function  on  an  element  A  vanishes  on  K  £  A,  then  either  K  -  A  or  K  is 

a  (one  dimensional)  curve. 

r 

(b)  Let  supp  4>p  ■  U  A^.  Assume  that  lA^I  >  lA^^  I .  k  -  l,...,r-l. 

k*l 

It  will  suffice  to  show  the  result  in  the  case  A'  *  A^. 

First  we  show  that  one  of  the  vertices  of  A^  must  be  P.  Suppose  that 
this  is  not  the  case,  and  let  Q  be  any  vertex  of  Ap  If  Q  is  proper, 
<|»p(Q)  ■  0.  So  at  least  one  vertex,  Q  say,  of  A^  must  be  improper  with 
<J»p(Q)  J4  0.  By  Lemma  3.1(b),  there  therefore  exists  an  element  A*  with  Q 
lying  inside  an  edge  of  A*  and  I  A*  I  >  |A^|.  Since  Q  €  A*  with  c|/p(Q)  J4 
0,  it  follows  that  A*  £  supp  (J^p.  But  this  contradicts  the  maximality  of 
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Next  notice  that  any  element  with  P  as  a  vertex  can  be  "connected"  to 
A^  either  directly  (If  (111)  of  Lemma  3.1(a)  applies)  or  Indirectly  by  way  of 
one  intermediate  element  (if  (11)  of  Lemma  3.1(a)  applies). 

We  now  prove  the  result  by  induction  on  k.  Suppose  that  we  can  "con¬ 
nect"  A^  to  Aj  for  all  k  <  j,  and  consider  If  P  is  a  vertex  of 

Aj+j,  then  by  what  was  said  above,  Aj^  can  be  "connected"  to  Aj.  If  P 
is  not  a  vertex  then  at  least  one  vertex,  Q  say,  of  A^  must  be  improper 
with  <j>p(Q)  +  0.  From  Lemma  3.1(b)  It  follows  that  there  exists  A*  f  V 
with  Q  $  A*,  1  A* l  >  IAj^jI  and  A*  (1  Aj^  an  edge  of  A.^.  Clearly 
A*  £  supp  4»p,  so  A*  ■  A^  for  some  k  <  j.  Thus  A.^  can  be  "connected" 
to  Aj  by  appending  A  to  the  connecting  chain  for  A^.  □ 


For  any  proper  node  P  we  shall  refer  to  Qp  ■  Interior  (supp(<pp))  as 
the  star  of  P.  We  shall  also  make  use  of  the  following  notation: 

*  {A  (  p:  A  5  supp(<J.p)), 

r*  *  U{y:  y  and  ®dge  of  some  A  €  Vp,  y  %  9Qp} , 

rp  -  aQp  n  a^o.i)2. 


The  three  tuple  <P*Dp»rp>  will  be  called  the  star  tuple  of  P.  Note  that 
the  star  uniquely  determines  Qp,  Fg,  rp  and  rp. 


Qp  -  Int(  U  a), 

P  A?Pp 


r*  -  U  (8A  -  Np). 

P  A€pp  P 


For  instance 


Furthermore  observe  that  whether  or  not  an  edge  contained  in  rp  is  primitive 
can  be  determined  from  Pp  alone. 

Let  us  also  introduce  the  notation 

F*  *  0{y:  y  and  edge  of  some  i  (  5,  y  %  3(0, 1)^}, 

r  -  r*  u  a^. 

Lemma  3.6. 

(a)  If  A  €  V  (and  V  contains  more  than  one  element),  then  A  €  Pp 
for  some  interior  proper  node  P. 

(b)  r*  -  U  r*. 

P 

(all  proper  nodes) 

(c)  blQ  -  u  ri. 

p  * 

(all  proper  nodes) 

Proof .  (a)  An  easy  induction  on  the  refinement  process  used  to  construct  the 

mesh  shows  that  for  any  element  at  least  one  vertex  is  an  interior  proper 
node. 

(b)  Clearly  ri  £  r*,  so  U  TS  2  T*.  For  the  reverse  Inclusion: 

P 

Suppose  y  is  an  edge  of  A  with  y  2  3Q.  Let  x  be  the  midpoint  of  y.  By 

Corollary  4(b),  £  <i»_(x)  “  1.  So  for  at  least  one  proper  node  P,  <J»_(x)  + 

P  F 

0.  It  follows  that  A  €  Pp  and  y  2  3Qp. 

(c)  Again  it  is  obvious  that  Ur*-  3*Q.  To  show  the  reverse.  Let 

P  P 

T  be  an  edge  of  A  €  p.  Then  by  (a)  of  the  lemma,  A  (  Pp  for  some  proper 
node  P.  a 

Suppose  y  is  a  straight  line  segment  which  when  extended  to  infinity  in 
either  direction  divides  R^  into  half  planes  and  x^  ^ .  Let  x  €  y 


and  suppose  v  *  (v^,v2)  is  a  sufficiently  smooth  function  defined  in  S+ 

*  x  n  {y:  |x-y|  <  e}  and  S“  -  x  fl  {y:  |x-y|  <  e}.  Let  n  and  n^  be 
the  unit  normals  to  y  pointing  out  of  x^+^  and  x^  ^  respectively.  We 

a 

shall  define  11  v* n ]]  at  x  by 


Ivnl 


lmt  v(y)»n^+^  +  lmt  v(y)«n^“^ 
y+x  y-*-x 

yes(+)  y€S^-^ 


Lemma  3.7.  Suppose  2  €  C^(Qp)  is  bilinear  on  each  &  €  Pp.  If 
[Vz»n]  *  0  on  each  y  €  rp  then  z  is  bilinear  in  Qp. 


Proof.  Let  A'  be  any  element  of  Pp.  The  function  z  is  bilinear  on  A'. 
Extend  this  bilinear  function  to  all  of  Qp,  and  call  this  extension  z*.  We 
claim  that  z  »  z*  on  all  elements  of  Vp.  To  show  this  it  suffices  by 
Theorem  3.5(b)  to  show  that  if  z  *  z*  on  A^  and  the  pair  Ap  A2  share  an 

edge,  y  say,  of  the  smaller  of  the  pair,  then  z  *  z*  on  A But  this 

follows  readily  since  (i)  z  is  bilinear  on  A2,  ( ii )  z  is  continuous 

across  y ,  and  (iii)  y  (  r£  and  so  the  directional  derivative  of  z  normal 

to  y  is  continuous  across  y.  a 

For  any  A  S  V  define  Q^,  the  influence  region  for  A  by 

■  Interior  (  U  supp  <J>p) . 

(A€0p) 

Introduce  the  notation  m  {a*  €  V’  A*  €  U  supp  4>p}  and  *  {y:  y 

A€Pp 

an  edge  of  some  A*  €  y  £  5^Q).  The  three  tuple  <A,P^,T^>  will  be 


called  the  Influence  tuple  of  A 


1.4.  We  shall  now  Introduce  a  restriction  on  the  "spread"  of  supp(<J>p). 

This  restriction  will  be  essential  for  our  analysis  later  on. 

For  any  set  S  S  R?  let  |S|  *  max  sup  (x. -y. |. 

1-1.2  x,y€S 

Let  K  >  1.  A  mesh  V  is  called  a  K-mesh  if  for  each  proper  node  P 

|supp(4>  )|  <  K  min  |&|.  (4.1 

MVp 

For  any  K  >  1  there  are  certainly  meshes  that  are  not  K-meshes.  For 
instance,  for  meshes  of  the  type  shown  in  Fig.  4(a),  supp(Cp)  is  always  the 


Figure  4(a).  Example  of  a  mesh  which  is  not  a  K-mesh. 

shaded  region,  yet  the  minimum  element  size  can  be  made  arbitrarily  small  by 
continuing  the  refinement  process  sufficiently  far.  The  K-mesh  property  is, 
in  some  sense,  only  a  local  property  as  it  permits  many  natural  forms  of  mesh 
grading.  For  instance,  the  grading  pattern  of  Figs.  3(b)  can  be  continued 
indefinitely  without  violating  (4.1)  with  K  *  2. 


Figure  4(b).  Example  of  a  K-mesh  with  refinement 


[It  Is  conjectured  that  the  above  definition  of  a  K-mesh  Is  equivalent 
to:  There  is  L  >  1  such  that  if  A, A'  €  D  and  A  PI  A'  *  0  then 

L_1  <  IA'|/IA1  <  L.  (4. 

By  equivalence  here  we  mean  that  If  (4.1)  holds  then  (4.2)  holds  with  L  ■ 
L(K);  and  If  (4.2)  holds  then  (4.1)  also  holds  with  K  ■  K(L).] 

In  what  follows  we  shall  always  assume  that  the  mesh  D  is  a  K-mesh. 
tie  now  state  some  Important  properties  of  K-meshes. 

Lemma  4.1.  Suppose  D  is  a  K-mesh,  then  there  exist  C  ■  C(K)  >  0  such 
that 

(1)  If  A  €  Dp,  CIAI  >  max  |A! . 

A€Pp 

(11)  If  P  is  a  proper  node,  card  Dp  <  C. 

(ill)  If  A  €  D,  card{P:  A  €  Dp}  <  C. 

(lv)  If  y  Is  an  edge  of  an  element  of  D,  card{P:  y  -  Tp)  <  C. 
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(v)  If  P  Is  a  proper  node,  card{Q:0p  H  0^  *  0}  <  C. 

(vi)  If  A*  €  0,  card{A:  A*  *  0  }  <  C. 

A 

(vli)  If  y  Is  an  edge  of  an  element  of  V ,  then  there  are  at  most  C 
nodes  on  y. 


Proof  * 


(1)  max  |  A I  <  Isupp  4»  I  <  K  min  IAI 
A€0p  ^  A£0p 


Isupp  (k,  I 

(ii)  card(Pp)  <  - y 

min  | A \ 


<  K2. 


(Ill)  If  A  €  0_,  then  Isupp  <|>_|  <  K  min  |A*I  <  K|A|.  So 

P  A*€0p 

U  (supp  4p)  £  where  Is  a  square  with  centre  at  the 

(A€Pp) 

center  of  A  and  side  2(R  -V2)  I A I  • 

Now  if  A'  €  U  0P,  then  A'  and  A  are  elements  of  the 
P  r 
(MVp) 

same  Vp  for  at  least  one  proper  node  P.  So 


I  A' I  >  j  Isupp  4«p  I  >  y  IAI 


lx* '  2  12 

Thus  card(  U  Vp)  <  — : - 7  <  4K  (K  -  y)  . 

P  (“IAI ) 

<A€0p)  K 

But  certainly  card(  U  V  )  >  y  card({  P:  A  €  0  })  since  P 

P  v 
(A 

must  be  the  vertex  of  at  least  one  element  of  Vp  and,  of  course, 
no  element  has  more  than  four  vertices.  Thus 


card{ Ps  A  €  0p)  <  16  K2(K  -  y)2. 
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(iv)  y  can  be  an  edge  of  at  most  two  elements.  By  (Hi)  each  of  these 
elements  is  contained  in  Pp  for  at  most  C  proper  nodes  P. 

Thus  y  can  be  contained  in  rp  for  at  most  2C  proper  nodes  P. 

(v)  By  (ii)  there  are  at  most  C  elements  in  Pp,  while,  by  (iii), 
each  such  element  in  turn  belongs  to  a  Pq  for  no  more  than  C 
proper  nodes  Q.  Thus  Pp  D  Pq  ^  0  for  at  most  C  proper 
nodes  Q. 

(vi)  By  (iii).  A*  is  contained  in  supp  4>p  for  at  most  C  proper 

nodes  P.  For  each  such  node  in  turn,  supp  4*p  can  contain  at 

2 

most  C  elements  A.  Thus,  card{A:  A*  €  P.)  <  C  . 

A 

2 

(vii)  If  y  is  on  9(0,1)  then  the  endpoints  of  y  are  the  only  nodes 

on  y.  Otherwise  suppose  y  is  an  edge  of  A.  Since  J  iL  ■  1, 

P 

at,  say,  the  midpoint  of  y»  then  for  at  least  one  proper  node  P, 


<|>p  >  0  at  the  midpoint  of  y»  In  fact  it  follows  that  <l»p  >  0 
everywhere  on  y,  except  perhaps  at  an  endpoint.  Any  node  on  y 
must  then  be  a  vertex  of  an  element  of  Pp.  By  (ii)  the  number  of 
such  elements  is  bounded.  □ 


We  would  expect  that  many  stars  Qp  are  essentially  identical  except  for 

a  translation  and  scaling.  Indeed,  as  we  shall  see  in  Theorem  3,  there  are  in 

fact  only  a  finite  number  (depending  on  K)  of  stars  up  to  translation  and 

rescaling.  To  set  the  scene  for  the  result  let  us  consider  an  arbitrary 

2  2 

proper  node  P.  Define  an  affine  transformation  xp:  R  -*•  R  , 


*p(x)  - 


[max  l A I 


(x  -  P). 


It  Is  readily  verified  Chat  xp(P)  -  0  and  max  |xp(A)l  •  1.  We  shall 

A€Op 

consider  xp  to  map  the  star  of  P  onto  some  "standard"  star.  The  mapping 
Tp  will  be  referred  to  as  the  star  transformation  of  P.  Let  us  set 


°0  “  W’ 


Q  -  { Tp(A):  A  €  I?p}, 


r;  - 


ro  -  r*  u  rj, 

*o  “  *p  0  tpX- 

The  two-tuple  <V$,  T^>  will  be  called  the  standard  star  tuple  of  P.  With  a 
harmless  abuse  of  notation  we  shall  refer  to  members  of  Vq  as  elements  and 
their  (closed)  sides  as  edges.  Notice  that  Qq,  rg,  Tq  and  Tq  can  be 
reconstructed  from  *n  lBuc^  t'ie  same  38  Qp>  r£,  rp  and  Fp 

can  be  expressed  in  terms  of  star  tuple  of  P.  Somewhat  less  obvious  is  the 
following  result. 


Lemma  4.2.  is  completely  determined  by 


Proof.  This  Is  clearly  the  same  as  saying  that  if  V  Is  a  mesh  and  P  is  a 
proper  node  of  V,  then  (|>p  is  completely  determined  by  <P,  p,rp>.  To  show 
this,  we  shall  construct  a  function  z  €  M(P)  using  only  our  knowledge  of 
<P»PD*r2>,  and  then  prove  z  ■  4s, : 


Arrange  the  elements  of  Pp  in  order  of  nonincreasing  size, 

say*  Suppose  that  z  has  been  defined  on  each  A;  for  j  < 

n.  To  define  z  on  An  it  suffices  to  specify  the  limiting  values  of  z 

from  within  A  at  each  vertex  Q  of  A  : 

n  n 

(I)  If  Q  -  P,  then  z(Q)  -  1. 

(II)  If  QM: 

(i)  If  Q  €  dQp,  then  z(Q)  -  0. 

(ii)  If  Q  e  Qp  and  Q  is  proper  node  (in  the  context  of  the 

<P,Pp,Tp>  this  just  means  that  whenever  Q  (  A*  for 
A*  €  Pp  then  Q  is  a  vertex  of  A*)  then  z(Q)  -  0. 

(iii)  If  Q  (  Qp  and  Q  is  improper  (just  as  in  (ii),  this  can  be  de¬ 

termined  from  Pp  alone),  then  by  Lemma  3.1(b)  Q  lies  in  the 
interior  of  an  edge  [P^^],  say,  of  some  unique  element  A* 

with  | A* I  >  lAnl.  In  our  present  case  A*  $  Pp  and  so  A*  *  Aj 

for  some  j  <  n  and  we  may  define  z(Q)  *  z*(Q)  where  z*(Q) 

is  the  limiting  value  of  z  at  Q  from  within  A*. 

This  enables  us  to  define  z  on  each  A  6  Pp.  Define  z  on  the 
remaining  elements  of  P  to  be  identically  zero.  By  the  above  construction 
z  has  a  well  defined  limiting  value  at  each  node  of  P,  and  so  by  Lemma 
3.2,  z€  ^(p)*  However  z  satisfies  (3.1),  and  so  by  Theorem  3.3,  z 
”  4>p-  a 

Theorem  4.3.  There  are  at  most  C  ■  C(K)  standard  star  tuples  <P q,Tq>. 

Proof .  It  will  suffice  to  show  that  there  are  at  most  C  possibilities 
for  Pq,  since  once  P q  Is  known  there  are  only  C  possibilities  for 
Iq.  (  Pq  has  at  most  C  elements  and  each  such  element  has  no  more  than 
four  edges.  Tq  must  be  the  union  of  some  of  these  edges.)  Since 


U  A I  <  K  min  |  A I  <  K, 


U  A  must  be  contained  In  a  square  with  centre  0  and  side  2K. 

Let  m  be  an  Integer  such  that  2  >  K.  Let  Q  be  a  square  with 
centre  0  and  side  2®*’*.  Define  a  uniform  grid  on  Q  of  grid  size  2_m. 

We  claim  that  each  A  s  Vq  consists  of  a  (whole)  number  of  grid  squares.  If 
we  can  show  this,  then  the  theorem  will  follow  at  once  by  simple  combinatoric 
considerations. 

Firstly  note  that  for  any  A€Pq,  certainly  |A|  ■  2  ^  for  some  Integer 
5  >  0.  Moreover, 

I A  |  >  min  |A'|  >  7I  U  A'l  >  2"“ 

A'€P0  *'<*><, 

and,  of  course, 

! A I  <  max  |A' I  ■  1. 

A'fP0 

In  particular,  If  (0,0)  Is  a  vertex  of  A,  then  A  must  consist  of  a  whole 
number  of  grid  squares.  To  prove  this  result  for  an  arbitrary  element  of  Vq, 
it  will  suffice  by  Theorem  3.5(b)  to  show  that  whenever  the  pair  A*,  A**  (  Vq 
are  connected  in  the  sense  of  Theorem  3.5(b)  and  A*  consists  of  a  whole 
number  of  grid  squares,  then  so  also  A**.  But  this  follows  readily  on 
recalling  (111)  of  Lemma  3.1.  □ 


We  also  shall  need  standard  forms  for  the  influence  tuples.  To  this  end 

2  2 

for  each  A  €  V  define  an  affine  transformation  t.:  R  ■+•  R 
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where  x.  Is  Che  centre  of  A.  It  Is  easily  seen  that  x.(A)  ”  [0,1]^.  Just 
A  A 

as  for  the  case  of  star  tuples  we  shall  consider  x^  as  mapping  the  influence 
tuple  of  A  onto  some  standards  reference  tuple.  Let  us  set 


°*  "  W* 

P*  -  {xa(A*):  A*  €  PA), 


The  two  tuple  <P*,TA>  will  be  called  the  standard  Influence  tuple  of  A ,  and 
clearly  Q*  can  be  reconstructed  from  the  standard  tuple: 

Q  -  Int(  U  A). 

A€P* 

Akin  to  Theorem  4.3  we  have 

Theorem  4.4.  There  are  at  most  C  -  C(K)  standard  influence  tuples  <P*,T*>. 

Proof .  It  will  suffice  to  show  that  there  are  at  most  C  possibilities 
for  p*,  since  once  P*  is  known  there  are  only  C  possibilities  for  T*. 
(By  (ii)  and  (iii)  of  Lemma  4.1,  P*  has  at  most  C  elements;  r*  must  be 
some  union  of  the  edges  of  these  elements.) 

To  prove  the  claim  for  p*,  we  apply  an  almost  identical  argument  to 
that  used  in  Theorem  4.3.  n 


1.5.  For  later  use  we  shall  need  to  know  something  about  the  approximation 
properties  of  the  finite  element  subspaces  flf(P)  and  Mq(P).  In  this  section  we 
shall  prove  some  results  in  this  direction.  The  methods  of  proof  are  more  or 
less  standard,  the  only  real  difference  arising  from  the  fact  that  the  star 


Qp  nay  spread  further  than  Just  those  elements  adjacent  to  P.  The  K-mesh 
assumption  however  restricts  this  spread  and  allows  us  to  retrieve  much  of  the 


standard  theory. 

Lemma  5.1.  Let  <PA,TA>  be  a  standard  influence  tuple.  There  is  a  constant 
C  -  >  0  such  that  if  z  *  Hr(Q#)  (r  *  1,2)  and  z  *  0  on  T* 

then 


inflz-ql  <  Clzl  (5.1) 

q€Q  r*Q*  r»°* 


where  Q  is  the  set  of  all  polynomials  of  degree  r  -  1  on  which  vanish 

on  IV.  By  1*1  (|»|  )  we  denote  the  usual  Sobolev  norm  (seminorm). 

r,w  r,u 

Proof.  Suppose  not,  then  we  can  find  z^  (n  ■  1,2,...) 


inflz<n)-ql 

q€Q 


n|z(n)| 


r,Q, 


(5.2) 


Without  loss  of  generality  we  may  as  well  suppose  that, 

inflz(n)-ql  -  1  (5.3) 

q€Q  r'U* 

(by  taking  suitable  multiples  of  the  original  z^  if  necessary)  and 

lz(n)l  <  2  (5.4) 

r»a* 

(by  adding  a  suitable  q11  €  Q  to  each  of  the  original  z^n^  if  necessary). 

From  (5.4)  we  can  conclude  using  Rellich's  Lemma  that  a  subsequence  of  the 

z(n)  converges  in  H^^Q^).  We  can  suppose  that  this  subsequence  is  the 

entire  sequence  (by  deleting  members  of  the  original  sequence  if  necessary). 

But  from  (5.2),  lim  |z^|  Q  *0.  Thus  z^n^  converges  in  Hr(Q*). 

Let  be  the  limit.  Obviously,  |z^"'l  Q  -  0,  and  as  q  is  connected 

r  t  w .  it 


(this  fact  follows  readily  from  Theorem  3.5(b)),  the  only  possibility  Is  for 

/a)  (cd  ) 

z  to  be  a  polynomial  of  degree  r  -  1  on  Q^.  Moreover,  z  ■  0  on  r* 

V  (os  ) 

(since  taking  traces  is  a  continuous  operation  in  H  (&*)).  Thus  z  €  Q. 
But  this  contradicts  the  limiting  form  of  (5.3).  a 


Theorem  5.2.  There  is  a  constant  C  >  0  such  that  for  any  z  (  Hr(Q)  f)  Hq 


(r  »  1,2)  and  any  K-mesh  V ,  there  exists  a  function  xz  €  such  that 


|z-*z|  <  CIAIr“slzl  , 

s,A  r,Q^ 


(s  -  0,1;  A  €  V ). 


(5.1 


Proof .  Let  p(t)  be  a  polynomial  in  one  variable  satisfying  /  p(t)dt  -  1 


and  /  tp(t)dt  *  0.  For  any  e  >  0  define 
0 


/  p(t/e) 

t  €  [0,e] 

♦£  ■  I  p(-t/e) 

t  €  ["6,0] 

(o 

Itl  >  e. 

P  of  V  set 

f  tain  |  A I 
z  ^ 

ME  p 


and  Sp  ■  {x^Q:  |x^  -  P^|  <  c?,  1  •  1,2}.  Our  choice  of  ep  guarantees 

1, 


Sp  s  Qp.  Por  any  y  (  HA(Q)  let 


yp  -  ^  y(*)#_  (*1-p1)*c  <*2”V  dx  /  4  dx 

P  ^  P  P 


(5.2 


and  define 
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Clearly  (Corollary  3.4(a))  icy  €  MQ(P). 

We  shall  now  prove  (5.1).  Suppose  A  $  V  and  let  <Z?*,rA>  be  the 
standard  Influence  tuple  of  A,  and,  as  usual,  write  Q*  ■  Interior  (  (J  A) , 


|ltyls,A  < 


l  IYpll<M 


P  P  s,A 


p^a°Q 


l  |Y  ||AI 
P 


1  a  cn  > 


by  Corollary  3.4  (b)  and  (c).  However,  from  (5.2) 


'V  «  <  cl4l'1|y'o,i 


whenever  A  €  Pp.  Thus 


|ltyls,A  <  C 


l  'y'o.Sp1^'8  <  c|A|’8|y,o,Q 


r 

1  Acf)  ) 


(5.3) 


by  (iii)  of  lemma  4.1  and  since  SpS  Qp  S  when  A  €  £>p.  Upon  rescaling 
to  the  unit  square,  (5.3)  becomes 


I  icy  °  x  .  I  o  <  Cly  or  | 

4.,(0,1>2  4  °’B* 


and  so 


|y  o  x  1  -  icy  o  t  1 1  <  Cly  o  t.1!  _  . 

A  A  s,(0, l)2  A  r’Q* 


(5.4) 
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If  C  is  a  polynomial  of  degree  r  -  1  on  which  vanishes  on 

I*.,  then  C  can  trivially  be  extended  to  all  of  Q.  Direct  calculation  shows 
that  *C  -  C*  Thus 


I*  °  “C.1  “  °  r  1 1  ,  -  inf  |  ( z+C )  o  T.1  -  x(z+C)  °  t.1!  , 

A  A  s,(0,l)2  C  *  *  8,(0, l)2 


-1. 


(5.5) 


<  inf I ( z+C )  0  t.  I 

r  A  r*Q* 


by  (5.4).  But  clearly  as  C  ranges  over  all  possibilities  C  °  ranges 

over  all  polynomials  of  degree  r  -  1  on  QA  which  vanish  on  r*.  Thus  by 
the  result  of  Lemma  5.1, 


infix  o 

C 


-1 


-COT 


-1. 


A  r,Qj 


<  C|  z  o 


TA  'r,Q . 


(5.6) 


Combining  (5.5)  and  (5.6),  and  rescaling  back  to  A  gives  the  desired  result 
(5.1)  upon  noting  that  the  constant  in  (5.6)  can  be  taken  to  be  mesh 
Independent  by  Theorem  4.4.  □ 


Corollary  5.3.  There  is  a  constant  C  >  0  such  that  for  any  K-mesh  V , 
(i)  If  z  €  Hr(Q)  fl  Mq  (r  -  1,2). 

Iz-xzI8>q  <  ChOO^Izl  ,  (s  -  1,2) 

(ii)  If  z  €  H^Q)  n  M0 

Proof. 

(i)  It  follows  directly  from  the  theorem  that 
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'*“**',  Q  <  ChCt?)1^8  (  I  Izl2  )/2 

*u  A  ’A 

<  ch(fl)r_ai*i  Q 

by  (vi)  of  Lemma  4.1. 

(ii)  If  A  €  Pp,  then  the  theorem  and  Corollary  3.4(b)  and  (c) 
shows 

(z-itz)H  <  C(  max  sup|D  <Jj  I II  z-irzBn  +  sup  |<J>  1 1  z-xzl  ). 
P  1,A  i-1,2  x€A  1  U,A  x€A  1,A 

On  the  other  hand,  If  A  ^  Op  then  obviously 

«ti)p(z-itz)«1  ^  -  0. 

Fixing  A,  squaring  and  summing  over  all  proper  nodes  P  gives 

l  lc|»_(z-itz))»?  <  CSzI?  n 

P  P  ,A  i,UA 

by  (ill)  of  Lemma  4.1.  The  result  now  follows,  just  as  in  (1),  from  (vi)  of 
Lemma  4.1.  □ 

We  shall  also  need  a  lower  bound  on  the  approximation  power  of  M(P). 
this  end  we  prove  the  following  lemma. 

Lemma  5.4.  Suppose  that  z  €  H*(Q)  and  that  there  is  an  open  disc  D  -  Q 
where 

(i)  z  €  C"(D), 

(ii)  max  inf  |D.  z(x)|  >  0.  (5. 

i-1,2  x€D 

Suppose  that  the  mesh  V  satisfies  for  some  k  >  0 


min  | A I 

A*D 

a£d 

Then  there  are  constants  C,  hg  such 

inf  lz-vl.  -  >  Ch(P). 

vdT(P)  1,0 

Proof .  Certainly 

>  J  inf  lz-wl? 

A€p  w  1  * 

(w  bilinear  on  A) 


>  £  inf  lz-wl2  . 

A€0  w  1 

A£D  (w  bilinear  on  A) 

On  each  element  A  £  D,  we  can  expand  z  as  a  Taylor  series  about  the 
centre,  x  say,  of  A, 

1  2  _ 

z(x)  -  w(x)  +  J  Dj^zCx)  (Xj-X^  +  r(x)  (x 

where  w  is  bilinear  on  A  and 

|r(x)|  <  C|x-x|3,  iD^r(x) |  <  C|x-x|2  (i  - 

for  some  constant  C,  independent  of  the  particular  A  under 
consideration.  Writing 


inf  Iz-vl,  _ 

xSW)  l'Q 


>  xh(P). 


that  if  h(P)  <  hr 


(5.8) 


z(x)  ■  w(x)  +  J  £  D  z(x)(x  -x  )2 
L  i-1 


we  obtain  using  the  finite  dimension  of  the  classes  of  functions  involved 


'l,A  »  C(j1IDll2U>l,l4|: 


w  bilinear  on  A 


>  C* I A  t  “ 


where  C*  >  0  is  independent  of  A  by  (5.7).  Hence 


Iz-wl  > 
1,A 


lz-wl  -  I r(x)i 

1  y  A  1  y  A 


(w  bilinear  on  A) 


(w  bilinear  on  A) 


>  C*  I A I  2  -  C  J  A 1 3 


>  CIAI' 


for  h(P)  <  hg,  and  hence  |AI ,  sufficiently  small.  Therefore  (5.9)  enables 
us  to  say 


inf  A  z- 


;-vlJ  >  C  l  IA! ‘ 


v€M(0)  A€P 


>  Ch(D)2  l  |A I2 


(5.10 


using  (5.8).  Note  however,  that  for  h(D)  small  enough,  D*  £  U  A  where 

A€P 


D*  is  a  disc  concentric  with  D  but  of  radius  half  that  of  D.  Thus 

I  IAI  >  area  D*,  and  the  lemma  follows  at  once  from  (5.10). 

A€0 


Let  us  remark  that  the  hypotheses  of  the  lemma  are  not  very  demanding  at 


2.  The  A-posteriorl  Error  Estimate 

2.1.  For  each  proper  node  P  of  p,  define 


Mp  ■  {v  (  H^(Qp):  v  -  0  on  Tp>. 


(Clearly  by  extending  functions  In  Mp  by  zero  to  the  remainder  of  Q  we  can 
consider  Mp  £  Hq )•  We  shall  associate  the  following  local  subproblem  with 
the  proper  node  P:  Find  r)?  g  Mp  such  that 


b(hp.v)  -  b(u-u,v). 


€  Mp. 


(1.1) 


After  an  Integration  by  parts  the  right  hand  side  of  (1.1)  can  be  written  out 
explicitly  In  terms  of  the  Input  data  of  (1.1.1): 


b(u-u,v)  -  l  j  (f-L(u))v  dx  -  J*  o(u)«n  v  ds 


Mflp  A 


+  /  (t-o(u)»n)v  ds. 

5'Q 


(1.2) 


0  2  ~ 

Note  also  that  as  long  as  PH  8  (0,1)  ,  then  €  Mq(P)»  and  that 

therefore 

b(u-u,^p)  -  0  (1.3) 

by  the  usual  finite  element  error  orthogonality  relation  (1.1.5b). 

We  are  now  able  to  state  and  prove  the  fundamental  error  estimate  of  this 
paper.  This  result  will  allow  us  to  estimate  the  global  finite  element  error 
lu-ul^  q  in  terms  of  the  solution  of  the  local  subproblems  (1.1).  The 
basic  ideas  behind  this  estimate  were  first  presented  in  [6]. 


Theorem  1»1.  There  is  a  constant  C  >  0  such  that  for  any  K-mesh 


lurul 


1,0 


< 


C(I  Ihpl 
p  r 


f/2 


(1.4) 


where  J  denotes  a  summation  over  all  proper  nodes  ?  of  V. 

P 

Proof,  (i)  The  right  hand  inequality:  By  the  coercivity  of  B  (see 
1.1.4c),  and  the  finite  element  error  orthogonality  relation  (1.1.6a),  we  have 
for  any  V  €  Mq(P) 


I  u-u  I 


2 

1.0 


< 

< 

< 


C  b( u-u, u-u) 

C  b(u-u,u-u-V) 

C  b(u-S,(][  <l>p)(u-U“V) ) 
P 


since  l  <|»_  ■  1  by  Corollary  1.3.4(b).  Thus, 

P 

lu-ul^  0  <  C  l  b( u-u , 4»p ( u-u-V ) )  . 

But  (|»p(u-u-V)  €  Mp  (certainly  4>p(u-u-V)  €  H^(Qp);  for  the  trace  behaviour 

2 

note  that  if  P  X  d(0,l)  then  <J>p  -  0  on  9Qp,  whereas  if  P  €  9Q  then 
u-u  -  V  ■  0  on  Tp  H  3^(0, l)^  while  d>p  ■  0  on  Tp  H  Q).  Thus 

b(  u-u,  4>p(  u-u-V))  -  b  ( rip,  (|;p(  u-u-V>) 


lu-ul^Q  <  C  |  |b(np,4>p(u-u-V))| 


‘  c  (I  'Via  j'2  )1/2 

P  ’  P  P  *  P 


and  we  have 


on  making  the  particular  choice  of  V  described  in  Corollary  1.5.3(11). 

The  left  hand  inequality:  Let  us  partition  the  set  of  proper  nodes 
of  V  into  disjoint  classes  Xj  (j  **  1,...,J(P))  say,  with  the  property 
that  whenever  P,P'  (  x^»  then  D  -  0.  It  follows  from  Lemma  1.4. l(v) 
that  we  may  always  ensure  that  J(V)  <  C,  for  some  mesh  Independent 
constant  C  >  0.  Now 

b(u-u,  l  t,  )  *  l  b(u-u,ri  ) 

P  P 


*  I  b(r)p,r)p) 


> 


c  l  'Vl.0p 


by  the  coercivity  of  B  over  Mp  c  Mq  (1.1.4c).  But 

J 


lb<u-u>  l  hp)l  <  Clu-ul .  Q  «  I  I  tlI,  c 

P  j.!  P  Xj  * 


2  M 


<  C  lu-uT  l  (  l  Jr,p»;  q  ) 

’  J-1  P  X  j  l,°* 

<  C  lu-ul  l  (I  IT,  i\  )V2 

*  j-1  P  ’  P 

<  C  l-u«1>Q  (J  IV?,q/2* 


The  left  hand  inequality  now  follows 


□ 
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The  significance  of  this  result  is  that  within  the  class  of  K-meshes  the 
ratio  6  “  (I  /lu-ul^  q  is  bounded  above  and  below  (away  from  zero), 

Independently  of  the  global  mesh  refinement  pattern  (as  long  as  the  K-mesh 
property  is  maintained).  Note  also  that  Theorem  1  demands  no  assumptions  on 
the  regularity  of  the  solution  u  (other  than  it  lie  in  Mg,  of  course). 

However  the  practicality  (1.4)  will  depend  on  whether  the  solutions  of 
the  local  subproblems  (1.1)  can  in  some  sense  be  effectively  estimated.  This 
is  the  matter  that  shall  concern  us  for  the  remainder  of  this  section.  Our 
main  result  will  be  Theorem  4.5  which  will  turn  out  to  estimate  the  Tip's  in 

^  A 

terms  of  the  "jump  residuals"  [o(u)«n]|  across  the  interelement 
boundaries. 


2.2.  In  this  and  the  following  section  we  establish  some  preliminary  results 
which  will  be  the  basis  of  our  effective  estimation  of  iTipl^  Q.  We  shall  be 
working  initially  with  a  version  of  the  subproblem  (1.1)  "standardized"  to  the 
standard  star  tuple  We  now  describe  this  standardized  problem:  Let 


Mg  -  {u  6  H1 (Qq) :  u  “  0  on  Tq} 


and  consider  a  continuous  bilinear  form  Ag!  Mg  x  Mg  -*•  R  which  satisfies  the 
conditions 


HAg(u,v) |  <  Plu«l>Qo»v,1>Qo 


-1...2 


I Ag(u,u) |  >  P  ,Ull,Qg 


u,v  €  Mg 


for  some  JJ  >  0.  Corresponding  to  any  Rq  6  I^CQq)  and  rg  €  ^(Fq)  let 
t)q  Mq  be  the  solution  of 

Aq(ti0.v)  -  /  RQv  dx  +  /  rQv  ds,  v  €  Mq.  (2.1 

Q0  ro 


Let  <Pq  €  Mq  be  a  function  with  the  properties 

(1)  4>0  >  0  on  Qq,  (2.: 

(11)  for  some  6  >  0,  <Pq  >  6  on  some  disc  contained  in  Qq. 

Finally  let  Gq  -  (Qq)  and  gQ  C  I^Cfg)  be  finite  dimensional  subspaces. 


Lemma  2.1.  There  is  a  constant  C  ■  C(Pq,Tq)  >  0  such  that  for  any  v  €  Mq 


IvIq  r  <  Clvlj 
u»1n  A»ur»» 


,0 


and  provided  Tq  t  0 


Ivl 


0,Qr 


<  CM 


l,Qr 


Proof .  The  first  estimate  is  a  standard  trace  result,  while  the  second  is  a 
Poincare-type  bound  on  the  L2  norm  in  terms  of  the  Hj  semi-norm.  a 


Lemma  2.2.  There  is  a  constant  C  ■  C(P,Pq,Fq)  >  0  such  that 


•Vi.Qq  <  c(,ro'o,Qq  +  ^o'o.rQ5, 


Proof . 


1 T10 1 1  ,Qn  <  C 


<  C  1/  RqTIq  dx  +  /  r0nQ  ds  I 


<  C(!RqI0  q  ItJqIq  q 


+  ,T1o,o,rnlrollo,rn^ 


By  Lemma  2.1,  ItiqIq  p  <  CItjqI  j  q  ,  and  of  course  "'no,0,Q  <  ,Tb'l  Q  ’ 

The  lemma  follows  Immediately.  □ 

Lemma  2.3.  Suppose  Gq  and  gQ  are  finite-dimensional  subspaces  of 
L2(Qq)  and  L2(Fq)  respectively,  then  there  is  a  constant  C 
-  C(0,Po,ro,Go,go)  >  0  such  that  if  Rq  €  Gq  and  rQ  $  g0  then 

*Vi,a0  >  0(,ro,o,q0  +  ''o'o.r,,)- 

Proof.  Suppose  the  result  does  not  hold,  then  there  are  sequences 

and  r^n\  (n  ■  1,...)  from  Gq  and  gQ  respectively,  and  corresponding 

solutions  of  (2.1)  such  that 

'"on>,i,an  <  k  (,Ro“>,o,a„  +  >•  <2-3 


We  may  without  loss  of  generality  assume  that 


,8o  'o>Qo  +  "ro  ' «0,r0  -  K  (2‘4 

In  particular  lR^n^Hg  ^  <1  and  Irg^lg  p  *  1»  By  the  finite  dimensions 

of  Gq  and  gQ  we  may  suppose  that 


Rj  ;  R  6  Gq  in  L2(Qq) 


♦  r  6  gg  in  ^CTq) 


Let  q  be 

the  solution 

of  (2.1)  corresponding  to  R  and  r. 

By  Lemma  2.2 

applied  to 

(n-r^0). 

<  C(»t— +  lr-r<")l0>1.0)  * 

0; 

but  by  (2.3),  tlQn)  ♦  0 

in  1*1.  n  and  so  we  conclude  that 

1»un* 

1 

o 

. 

Thus  for  any  v  €  Mq 

0  -  AflCrjjV)  "  /  Rv  d x  +  /  rvds. 

« 


But  this  can  only  mean  that  R  ■  0  and  r  ■  0.  However  this  clearly 
contradicts  the  limiting  case  of  (2.4).  a 

Lemma  2.4.  Let  Gq  and  gg  be  as  in  Lemma  2.3,  then  there  is  a  constant  C 
«  CCQqjTqjGq ,gg ,<pQ)  such  that  if  Rq  €  G0,  r0  ?  go  and  p  €  PQ  (the  set  of 
constant  functions  on  Qq) 

'Vo.q  <  c(lrolo,r0  +  IRo“plo,o  +  *0  0  dx  +  /  ro  o  d8|)* 

U  U  U  Uq  Tq 


Proof.  We  shall  again  argue  by  contradiction.  So  suppose  the  result  is  not 
true,  then  there  are  sequences  R^n\  and  p^n^  from  Gq,  gQ  and 

Pq  respectively,  such  that 


IR, 


>  n(lr£n)l  +  IR^n)  - 

°,rn 


0,Qr 


p(n)l 


0,Qr 


1/ 

Q 


<n\ dx  +  / 

0  ro 


.(»).. 


ds  I ) 


(2.5 


where  we  can  suppose  that 


By  the  finite  dimension  of  Gq  we  may  as  well  assume  that  Rq  R  €  Gg  in 

L2(Cq).  But  by  (2.5),  IR$n)  -  p(tl)«0ffl  *  °*  *»«*■  p(n)  ♦  R  in  L2(QQ). 

But  Pq  is  certainly  a  closed  subspace  of  L2(Qq)  (since  it  is  finite 
dimensional),  and  so  R  €  Pq,  l.e.,  R  is  a  constant  on  Qq.  Now  (2.5)  also 
shows  that  *  0  in  I^^q)  and  that 

I  /  R^n^tpQ  dx  +  /  rQa^Q  dsl  >  0  as  n  -*■  °».  Thus  we  must  have 

flo  ro 

0  ■  lim  /  r£QVq  dx  ■  /  R<Pq  dx. 

n-H»  Qq 

The  only  way  this  can  happen  when  <Pq  has  the  properties  (2.2)  and  R  is  a 
constant,  is  for  R  **  0.  But  this  contradicts  the  limiting  case  of  (2.6).  c 

Lemma  2.5.  There  is  a  constant  C  *  C(Qq,Pq)  >  0  such  that  for  any  z  €  8Q  “ 
(z  €  C°(Qq):  z  bilinear  on  each  A  €  Pq}, 

inf  IDz-pln  _  <  C{(  X  |Dz  |?  )1/2  +  (/  |Vz«n  B  daf'2  }  (2.7) 

°  ’  T°0 

where  Dz  denotes  any  (element-by-element )  zeroth,  first  or  second  order 
derivative  of  z. 


Proof.  Let  M  «  {z  (8  :  Dz  €  Pq}*  Clearly  the  quotient  space  Bq/M  is 

finite  dimensional,  and  (2.7)  will  follow  if  we  can  show  that  both  sides  of 

(2.7)  define  norms  on  8q/m*  In  both  cases  it  clearly  suffices  to  show  that 

if  the  respective  quantity  in  (2.7)  vanishes  then  z  €  M: 

(i)  The  left  side:  If  inf  IDz-p!n  n  -  0,  then  Dz  is  a  constant  on 

PSP0 

Qq  since  Pq  is  finite  dimensional  and  hence  closed. 


(ii)  The  right  side:  If  (J  |Dz|,  )  ^  «  0,  then  Dz  is  constant 

*»A 


on 


A€R 


each  A££L.  On  the  other  hand,  if  J  HVz»hu  As  ■  0,  then  tvz*nll 


0  on  Tq,  and  so  by  Lemma  1.3.7  z  is  bilinear  on  Q^.  In 


particular,  Dz  €  Cu(Qft).  The  net  result  being  that  Dz  €  Pq. 


We  shall  now  relax  the  requirements  of  Lemmas  2.3  and  2.4  that  Rq  and 
rQ  lie  in  the  finite  dimensional  spaces  Gq  and  gQ.  Instead  we  shall 
require  that  they  can  be  approximated  from  within  these  spaces. 


Lemma  2.6.  Let  Gq  and  gQ  be  as  in  Lemma  2.3,  and  suppose  that 


inf  IVR*«n  o  +  lnf  lrn-r*«ft  r  <  e. 
R*€Gq  v  *U0  r *(gQ  U,r0 


Then 


(i)  there  is  a  constant  C  -  C(0,QQ,rQ,GQ,gg)  >  0  such  that 


c  ('Vo,q0  +  ,ro"o,r0"e^  <  ,Tb "l ,cQ  <  c^,Ro“o,c0  +  nro»o,r0^ 


(2.1 


(ii)  there  is  a  constant  C  ■  C(Qq  ,Vq  ,Tq  ,Gq  ,gQ  ,<|;q)  >  0  such  that  if 


/  Vo  dx  +  /  Vo  ds  *  0 


(2.1 


then 


IR, 


■0'0,Q0  <  C^r°,0»r0  +  p“p  "Vp,0,Q0+e^ 


(2.1( 


Proof.  Choose  R*  €  Gq  and  r*  6  g0  such  that 
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IRn-R*ln  n  +  lrn-r*l!n  „  <  2e.  Let  h*  be  the  solution  of  (2.1) 

u  U»Q0  u,10 

corresponding  to  r*  and  R*.  Applying  Lemma  2.2  to  the  difference  (tu— p*) 
we  get 


<  «"o-4"o,q(|  *  ,ro-r*'o,r0 


)  <  Ce, 


Using  Lemma  2.3, 


"Vo.Gq  +  "ro"o,rc 


<  KR*»n  +  I r*ln  _  +  2e 
u»“n  u*1n 


<  CIr)*i,  +  2e 


<  C(lr)0llfQ  +  e), 


proving  the  left  side  of  (2.8).  The  right  side  follows  from  Lemma  2.2. 


For  (2.10)  we  have 


'Vo,o0  <  <*"o,a0  +  2‘ 


while  Lemma  2.4  gives  for  any  p  €  Pi 


lQ  Q  <  C(  Hr*B0  +  lR*-pfl0  +  1/  R**q  dx  +  /  r*<l0  ds  I ) 

,U0  u,10  ,U0  Qn  Tn 


<  c(,ro!o,rn  +  ,Ro”p*o,an  +  ^  Vo  dx  +  [  Vo  dsl  + 

U  0  a0  ro 

<  c(,ro"o,r0  +  "Ro'P"o,Q0  +  E) 


by  virtue  of  the  orthogonality  relation  (2.9) 


□ 
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2.3.  In  Che  last  section  we  considered  a  "standardized"  form  of  the  local 
subproblems  (1.1).  In  this  section  we  extend  those  results  to  the  actual  star 


tuple  <p,  p,rp>. 

For  any  proper  node  P  with  standard  star  tuple  <t?Q,rg>,  let 
2  2 

x  :  R  be  the  star  transformation  of  P  (see  §1.4),  that  is 

P 


x  (x)  -  ±  (x-P) 

p  p 


where  p  *  p_  *  max  |A|.  (Notice  that  p  <  1.) 

P  A*Pp 


Corresponding  to  the  finite  dimensional  spaces  Gq  and  gg  of  §2.2  let  us 
define 


Gp  -  {v  $  L2(Qp):  vo  €  GQ}  c  L2(Qp), 


gp  -  (v  ^  L2(rp):  v  e  Xp1  €  Xp^g}  c  L2(rp). 


At  this  stage  it  is  convenient  to  consider  a  more  general  version  of  (1.1): 
Find  t)  €  Mp  such  that 


b(r|,v)  *  /  Rv  dx  +  /  rv  ds,  V  v  €  Mp  (3.1) 


where  R  €  L^Qp)  and  r  €  L2(Tp). 


-v 


Lemma  3.1. 


c  ■  Pp  inf  1R-Rln  _  +  Pp2  inf  «i 


•'o,rp  <3-2> 


%  /.  %  *.  /.  •-  ;.  •• 


(i)  there  is  a  constant  C  -  C(P0,r0,G0,g0)  such  that 


c‘lt<’p,R,o,Qp  +  pp2 lrlo,r."s)  ‘  <  c(Pp*"o.o.  +  pp2 lr"o,r J 


(ii)  there  is  a  constant  C  -  C(t?0,rQ,G0,g0)  such  that  if 


/  R4«_  dx  +  /  r<J)_  ds  -  0 

Q  r 

WP  *P 


Pp,R,o  Q  <  c(pp2  ,r,o  r  +PP  inf  ,R-P*o  Q  +£) 

,BP  p  0,rp  p  P€P0(Qp)  ’  p 

where  PqCQj)  is  the  space  of  constant  functions  on  Qp. 

Proof.  It  is  clear  that  Mg  »  {v  f  H * (Oq ) :  v  o  i  €  Mp}.  Define  a  bilinear 
form  Aq:  Mg  x  Mg  +  R  by 


Aq(u,v)  -  b(u  o  Tp,  V  O  Tp), 


It  follows  from  the  boundedness  of  b(»,»)  on  Mp  x  Mp 


|An(u,v)|  -  lb(u  o  t  vo  tv)I 


<  ciu  o  T_| ,  _  Jv  o  T 


P1  l,Qplv  °  Vl.Qj 


C(p2,U,5  Q  +  |U|?  Q  )l/2(p2«V»0  Q  +  |Vll  Q  )l/2 

u,wo  ,uo  u,uo  ,ao 


<  c lul.  Q  Ivlj  Q  . 
i,uQ  i,Uq 


In  addition,  by  the  coercivity  of  b(»,»)  over  Mp, 


•  -  •*  »  '*  -  *  ,*  ,  *  .  *  «•  ,  «  •  *  .  *  *  '  ■»  *  v  ‘  *«*-•-•  .  •  .  *  .  .  •* 


|Ag(u,u)| 


|b(u  o  Xp,  u  o  Tp) |  >  alu  o  Tp»^  £ 


>  C|u| 


l.Or 


>  Ciul 


l.Qf 


using  lemma  2.1,  noting  that  Tq  #  It  is  readily  verified  by  a  change  of 
variable  of  integration  that  for  any  v  f  Mg 


AqOh  °  *piv)  -  b(ri,v  o  t"1)  -  /  R(v  o  t"1)  dx  +  JT  r(v  o  r'^ds 


-1  ■ 


-I- 


/  (p  R  o  T  )v  dx  +  /  (pr  O  Tp ) v  ds. 


So  we  are  in  a  situation  covered  by  the  analysis  of  §2.2  with  t7q  *  rj  °  rp. 


lg  -  p  R  o  Tp  and  rg  -  pr  o  Tp.  Furthermore,  for  any  S  f  Gp  and  r  €  gp, 


2  s 


we  have  p  R  o  Tp  f  Gfl  and  pr  o  Tp  (  gQ  with 


»p2R  o  Tp  -  p2  R  o  TpA0>Qo  +  ipr  o  t?  -  p'r  o  Vo.r 


0 


p«**«  +  pl/z  ffr”^no,rp' 


By  virtue  of  (3.2)  we  are  now  able  to  apply  Lemma  2.6  to  infer  from  (2.8) 


C^'P2ROTP'o)Bg  +  lproTp«0)r()-e)  <  "^Vl.Qg 


(3. 


<  C(lp2R  o  ^p»o>Qo  +  "Pr  o  Tplg^g)* 


Upon  rescaling  back  to  the  actual  star  Q_, 


,T1°Vl,Qn  "  ,T1  °  TP!0,Qn+  |T1  °  XPll,Q0  "  -2  |Tll0,Qp  +  ,T>' 


>  Itllt  n  J 


while  by  Lemma  2.1  we  also  have  the  opposing  bound 


‘  °">°Vi,a0  ‘  c"»i,op  *  c""i,ap- 

Rescaling  the  other  terms  in  (3.3)  leads  at  once  to  (i)  of  the  lemma. 

To  show  (ii)  note  that  ■  4>p  o  xp  satisfies  (2.2)  and  that 

/  (p2R  o  +  /  (pr  o  xp)(J>0  -  /  R<t>p  dx  +  /  r(j,p  ds  *  0. 


Part  (ii)  of  Lemma  2.6  then  gives 


«p2R  o  x  I  <  C(«pr  o  x  B  +  inf  Bp2R  o  x -  PB  +  e), 

p  0,Qo  P  U»F0  P€Pq  P  U»Q0 

where,  since  is  determined  by  the  standard  tuple  <Qq,Pq,Fq>,  the 

constant  C  of  Lemma  2.6  (ii)  is  in  fact  C  *  C(QQ,0Q,rQ,Gg,gg) .  Rescaling 
to  Qp  gives  the  result.  □ 


2  2 

2.4.  We  now  return  to  the  task  of  effectively  estimating  (£  BripB,  Q  )  .  Th 

P  P  i,UP 

local  subproblems  (1.1)  fit  into  the  pattern  of  §2.3  if  we  set  (see  (1.2)) 


R  *  (f-L(u)) 


on  each  A 


-|Ia(u)»n]]  on  r* 


»  _  ,7\  -  - 


on  a^O,!)2. 


To  complete  the  framework  required  for  applying  the  results  of  §2.3  we  define 


Gq  *  {v  €  L^CQq):  v  is  constant  on  each  A  €  Pq}, 

gg  *  {v  €  I^CTq):  v  is  linear  on  each  primitive  edge  y  Q  Tq} • 

It  is  clear  that  Gq  and  gQ  are  finite  dimensional  and  that  they  are 
determined  completely  by  the  standard  star  tuple 


Lemma  4.1; 


(i) 

l  IAl20ullf  -  o(l) 

A*p 

as  h(P ) 

-*•  0. 

(ii) 

I  Pp  inf 

Sr-RIq  - 

o(h(P)2) 

as 

h(P)  *  0. 

P  *  R€Gp 

(iii) 

I  PP  inf 

«r-?ll2  - 

o(h(P)2) 

as 

h(P)  0. 

*  *€gp 

P 

(iv) 

There  exists  r*  €  L2(T)  which  is 

linear 

■  on  each  primitive  edge  y  £  T, 

which  agrees  with  r  at  all  nodes 

lying 

on  T,  and  which 

satisfies 

I  IAI 
A 

#r-r*B  “ 

aAnr 

o(h(P)2) 

as 

h(P)  -►  0. 

Proof.  To 

prove  (i) 

t 

T  !A|2»uJ2  .  < 

h2(P)llul 2 

+  l 

IA|2|u|2  . 

A€Z7 

1.0 

A€P 

By  the  projection  property  (1.1.5)  of  the  finite  element  solution  the  first 
terms  obviously  tends  to  0.  For  the  second  term 


I  IA| 2 |u| 2  A  <  C  l 
A€J?  A€p 


|u-v|! 


(v  linear  on  A) 


9  9 

<  C  £  (lu-ul.  +  inf  lu-vi.  .) 

A€p  *A  v  l,A 

(v  linear  on  A) 


<  c(lu-ul^  +  £  inf  lu-vll?  .) 

1,Q  A  v 

(v  linear  on  A) 


which  certainly  tends  to  zero  as  h(P)  -*•  0. 

Now  turn  to  (ii)  of  the  lemma.  On  each  element  A  €  flp,  R  is  a  sum  of 
terms  each  of  the  form  X.(x)£(x)  where  k  has  bounded  derivatives  of  all 
orders  on  Q  and  £  -  l,u,D^u  or  D^u  (i,j  “  1,2).  We  can  find  a  constant 
function  0  on  A  such  that 


.X5-ei0>A  <  ciai«xc«1>a 


<  CIAHCI. 


<  C|A|(IuI2^a  +  I A I )  ■ 


inf  IR-RI.  -  C{  l  OAI2^  .+  IA|4)f/2, 


and  so 


I  p2  inf  #R-RIq  <  C  l  p2  l  (|A|29uI2  A  +  IAI4) 
P  r  R€Gp  ,UP  P  A€Pp 


<  C  h(0)2(  l  | A 1 2HuI?  .  +  h(P)2) 


where  we  have  made  use  of  the  fact  that  each  element  A  €  P  belongs  to  a 
Pp  for  at  least  one,  but  no  more  than  C  -  C(K),  proper  nodes  P  (Lemmas 
1.3.6  and  1.4.1).  Part  (11)  of  the  lemma  now  follows  Immediately  on  using 
(1). 

For  part  (111),  observe  that  on  each  element  edge  y  making  up  rp,  r 
Is  a  sum  of  terms  of  the  form  X(t)5(t)  where  \  has  bounded  derivatives  of 
all  orders  (independent  of  P)  and  5  -  1  or  a  limiting  value  of  D^u  (1  * 
1,2)  on  the  edge.  There  is  no  loss  of  generality  in  assuming  tha  D^u  is 
linear  on  this  edge.  We  can  find  a  linear  function  p  on  y  satisfying 

l«-pi0  v  <  C  | y  | 2 1 -  v 

<  ciyi2(»5I2>y) 

<  ClY|2(  l  ID  ul  +  |y|1/2) 

1-1  1  i,Y 

where  D^u  Is  the  limiting  value  on  y  from  within  an  element  A  €  Pp.  We 
have 

ID  .ul .  <  ClAl"  1/2  lul-  . 

i  1,Y  2, A 

and  so 

l\5“PI0  <  C|A|1/2  (|A|IuS2jA  +  |A|2).  (4.2) 

Therefore 

inf  lr-rlQ  <  Cp^2  {  J  |A|2lu»2  +  |A|4}1/2 

r€gp  U’rP  P  A<tfp  Z’A 

and  so,  much  as  part  (1),  squaring  and  summing  over  all  proper  nodes 
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i 


l  Pp  inf  ,r-f,Z0  r  <  C  l  Pp  l  ( IA I 2«ul ?  .  +  IAI4) 
P  r€g„  U,rP  P  *  A€P„ 


P  r€g. 


<  Ch(0)*Q  IA|2lul2  +  h(D)Z). 

A 


Again,  on  using  (i),  part  (iii)  of  Che  lemma  follows. 

Parc  (iv)  of  Che  lemma  likewise  follows  from  (4.2)  and  (4.1),  since  we 
may  as  well  choose  0  Co  inCerpolaCe  X.£  ac  Che  endpolnCs  of  \.  □ 


We  are  now  able  Co  sCaCe  and  prove  Che  main  resulCs  of  Chis  paper 
concerning  an  equlvalenc  esCimaCor  for  Che  energy  norm  of  Che  flnlCe  element 
error.  Because  of  Che  coercivlcy  of  b(*,*)  there  is  no  need  Co  distinguish 
between  equlvalenc  estimators  for  (b(u-u,u-u) fa  ,  0r  lu-ul^  For 

definiteness  we  shall  phrase  our  result  1  »  terms  r>f  the  Latter. 


Theorem  4.2.  There  is  a  constant  C  >  0  such  that  for  any  K-mesh 


lu-ulj  q  <  CE 


where 


(4.4a) 


E  -  (  L[HI2«m5,4+  '*"<34nr)1/2 


If  in  addition  for  some  L  >  0 


lu-ul^  g  >  Lh(D) 


(4.5) 


then  there  are  constants  C  >  0,  hg  >  0  such  that  for  any  K-mesh  V  with 

h(0)  <  hQ 


C  aE  <  lu-ul^  g  . 


(4.4b) 


*_  -  a.*  .*  «.*  t  ^ ' *  a_  •  *  a*  a.  *  ^  t*  «'  •'  •*  •_  «  -  »-  .  -  ■*- 


By  Lemma  3.1 


T roof .  Consider  a  particular  proper  node  P  and  its  star  Qp. 

C-1(E  -p  inf  IR-RIq  +  Pp2  inf  Ir-rt  )  <  It,  I,  <  C  Ep  (4.6) 
P  R6Gp  °*°P  P  r€gp  °»rP  P  i,QP 


where 


l 

A€Pp 


lAUrl 


2 

0,6A 


Here  we  have  used  the  fact  that  for  any  A  6  0p,  C~^pp  <  IAI  <  pp  (see  Lemma 
1.4.1). 

The  constant  in  (4.6)  depends  only  on  the  standard  star  tuple 
associated  with  the  proper  node  P,  but  is  otherwise  independent  of  the  mesh 
or  the  node  P.  However,  by  Theorem  1.4.3,  for  any  K  mesh  there  are  at 
most  C  ■  C(K)  <  <®  possible  standard  star  tuples.  We  may  therefore  regard 
the  constant  C  in  (4.6)  as  uniform  for  all  proper  nodes  P  and  all  meshes 
of  the  specified  class. 

i  o  o  o 

We  claim  that  C-1  EZ  <  £ep  <  C  E  •  To  see  this  we  need  only  recall 

P 

(i)  any  element  A  (  V  belongs  to  a  Qp  for  at  least  one  but  no  more 
than  C  •  C(K)  proper  nodes  P  (see  Lemmas  1.3.6  and  1.4.1) 

(ii)  any  edge  contained  in  T  is  contained  in  a  Tp  for  at  least  one, 
but  no  more  than  C  *  C(K)  proper  nodes  P  (see  Lemma  1.4.1). 

Thus  squaring  (4.6),  summing  over  all  proper  nodes  P,  and  using  (ii)  and 
(ill)  of  lemma  4.1  gives 


C_1  (E2"h(D)2o(  1 ) )  <  l  ItuI?  0  <  C  E2.  (4.7) 

p  r  ‘.Up 


Theorem  1.1  allows  us  to  conclude  from  (4.7) 


Theorem  4.2  requires  no  extra  assumptions  on  the  regularity  of  the  exact 
solution  u  other  than  (4.5).  However  this  requirement  is  hardly  a 
restriction  at  all  being  satisfied  in  all  but  trivial  cases  (see  Lemma  1.5.4). 

Although  the  estimate  of  Theorem  4.2  is  sharp  in  the  sense  described 
above,  it  can  be  simplified  even  further.  It  turns  out  that  the  L2  residual 
contribution  can  be  dropped  from  E,  giving  an  expression  solely  in  terms  of 
"line  residuals,"  while  still  maintaining  the  two-sided  estimate  (4.4).  This 
will  be  shown  in  Theorem  4.5.  It  is  natural  to  ask  whether  the  "reverse" 
simplification  is  valid,  that  is,  can  the  "line  residuals"  be  dropped  from  e 
without  affecting  the  estimate.  It  is  not  difficult  to  see  that  this  cannot 
be  true.  Consider  the  simple  case  where  L(«)  is  the  Laplaclan.  Then  for 
bilinear  elements,  L(u)  ■  0,  and  the  L2  residual  becomes  just 

(I  1 A 1 2  /  If  1 2  dxf/2  >  Ch(0)lfl 
A  A 

if  the  mesh  V  is  quasiuniform.  However,  if  the  solution  u  has  some 

o  _i* 

singular  behaviour  so  that  u  if  H  (Q)  then  lu-uJj  ^  must  converge  at  a 
slower  rate  than  0(h(D)).  So  the  L2  residual  alone  cannot  suffice  for  E 
in  (4.5)  in  this  case. 

For  proving  Theorem  5  we  shall  need  the  following  lemmas: 

Lemma  4.3.  There  is  a  mesh  independent  constant  C  >  0  such  that  for  any 
proper  node  P  and  any  z  c  8(Qp)  ■  (z  €  C^(0p):  z  bilinear  on  each  A^0p}, 

inf  IDz-plQ  <  C{p  (  l  I Dz l2  /2  +  Pp(/  [[Vz»n  B2  ds/2} 
p€P0(Qp)  °’°P  ?  A?Dp  p  r* 

where  Dz  denotes  any  (element-by-element)  zeroth  (s  *  3/2),  first 
(s  -  V2  )  or  second  (s  -  -  V2  )  derivative  of  z. 


Proof.  This  is  just  the  rescaled  version  of  Lemma  2.5,  taking  note  again  that 


there  are  only  a  finite  number  of  standard  star  tuples.  □ 

Lemma  4.4.  Suppose  y  €  Q  is  a  straight  horizontal  or  vertical  line  segment, 
then  for  any  z  €  $(P)» 

/  |[Vz*n  B2  ds  <  C  /  |o(z)»n)  J  2  ds. 

Y  Y 


Proof.  For  definiteness  suppose 
that  y  is  vertical  and  let  super¬ 
script  (-)  and  (+)  indicate  limits 
from  the  left  and  righthand  sides 
of  y  respectively.  (See  Fig.  5.) 

Since  all  the  a^j^*  are  continuous 

in  Q,  Figure  5.  The  scheme  of  notation. 

I<r(z).;i  -  l  (<»1:jk;,Dtzjnk)<+)  +  (•1Jkj,°1*Jnk)<’>t 


I  {aiju['(Vj)(*  +  11 


i.j,* 


since  n£  *  0,  n[+^  ■  -1  and  n[~^  *  1.  But  * 8  continuous  across  y, 


ff<r(z)»nB  -  l  a....  [-(D,z,r+)  +  (D  z,)(_)] 


But  the  matrix  a11  is  uniformly  positive  definite  on  Q  (see  1.1.4b)  and 
so  det(a^)  >  C  >  0  on  Q.  It  follows  then  that  on  y 

||[Vz«n]]|  <  C|  llo(z)  «n]|  | , 

and  the  lemma  is  proven  in  this  case. 

The  case  of  y  horizontal  is  treated  similarly.  □ 

Theorem  4.5.  Suppose  that  for  some  L  >  0 

lu-G^  q  >  Lh(P) 

then  there  are  constants  C,hg  >  0  such  that  for  any  K-mesh  V  with 
h<0)  <  hg 

C_1E  <  lu-ul  <  CE  (4.! 

1.0 

where 

E  ’  I*  l4"r,o,84nr^-  <4-' 

Proof.  First  note  that  the  left  hand  side  of  (4.9)  follows  trivially  from 
(4.4b)  since  E  <  E.  To  prove  the  right  hand  side  it  will  suffice  to  show 

l  |A|2IRIq  <  E2  +  o(h(D)2)  (4. 

A  * 

with  the  o(l)  term  valid  as  h(P)  +  •. 

For  any  Interior  proper  node  P  we  have  (see  (1.2),  (1.3)) 

/  R<|)  dx  +  J  ri  ds  ■  0 

A  *  TV  * 


with  R  and  r  as  in  (4.1).  Thus  we  may  apply  part  (ii)  of  Lemma  3.1, 


Pp'R,0,Qt 


<  C(pl^  lrln  _  +  p„  inf  IR-plft  +e_),  (4.11) 


,p  —  0r  kp 


?  wVV 


0,Qb  tP> 


where  Pg(Qp)  is  Che  space  of  constant  functions  on  Qp  and 


e  -  p  inf  »R-RIQ  +  p^inf  lr-SlQ  _  . 
R*G„  ,aP  resu  ,XP 


We  shall  concentrate  on  the  term  inf  IR-pl n  n  for  a  while.  As 

peP0(Qp)  u’qp 

the  coefficients  of  L(»)  are  smooth  on  Q,  let  us  write  £(•)  for  L(») 
but  with  its  coefficients  replaced  by  their  averages  over  Qp.  If  we  then  set 

R  -  l  (F-L(u)), 

we  have 


inf  IR-pl n  <  !R— Rln  +  inf  IR-pln 

p€P0(Qp)  °*Q?  °’QP  P€PQ(Qp)  °’QP 


<  Cp  lul,  +  inf  lR“pl0  q 
‘  2’°P  p€P0(Qp)  U,CP 


(4.12) 


<  Cp  (lul  +  IQ  I)  +  inf  I  l  L(u)  -  pl0 
2>aP  P  p€P0(Qp)  A€Pp 


where  lu^  q  is  to  be  understood  in  an  element-by-element  sense. 

’  P 

Now  L(u)  is  just  a  sum  of  terms  of  the  form  X(Du)  where  X  is  a  constant 
and  Du  denotes  a  zeroth,  first  or  second  derivative  of  the  funtion 
u  €  8 (Qp ) •  (The  constant  X  is  bounded  independently  of  P  and  0).  We  may 
clearly  apply  Lemma  4.3  to  obtain 
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inf  I  l  L(u)  -  plQ 
>«(0  AeD„ 


0  P  p 


C(p plul_  +  p  ^2  (/  |7u»n  ^ds)1^2) 


P  2  Q  Kp 

r  *»“p  r  P* 


(4.13) 


<  C(p_lul _  +  p~  ^  (/  |o(u)»n  l2  ds)1/2 


P  2  Q  KP  VJ 

,up  ?  r* 


by  Lemma  4.4. 


Substituting  (4.13)  in  (4.12),  and  then  this  in  turn  into  (4.11)  gives 


PP,R,0  <  «p^.r.  +«*) 


(4.14) 


where 


<  ep  +  Pp(lul2>Qp  +  IQpl). 


Squaring  (4.14)  and  summing  over  all  interior  proper  nodes  P  gives 


I  |A|2IRI2  <  C  I 

A€D  ,A  P 


PP,R"0,QT 


(P  an  interior 
proper  node) 


(4.15) 


<  C(E2  +  l  (e_)2  +  h(D)2(£  |A|2lul,  +  h(D)2)) 
P  r  A 


where,  much  as  in  the  proof  of  Theorem  4.2,  we  have  used  the  following 
properties  of  our  meshes: 

(i)  for  A  €  Dp,  C  *Pp  <  |A|  <  Cpp  (see  Lemma  1.4.1) 

(ii)  any  element  A  (  V  belongs  to  Dp  for  at  least  one  interior 

proper  node  P  (see  Lemma  1.3.6)  and  at  most  C  *  C(K)  proper 
nodes  P  (Lemma  1.4.1) 


>•„»  „»  v ' 


(iii)  any  edge  contained  in  y  is  contained  in  Tp  for  one  more  than  C 
*  C(K)  proper  nodes  P  (Lemma  1.4.1).  □ 

The  desired  result  (4.10)  follows  at  once  from  (4.15)  and  the  assumed 
lower  bound  on  lu-ul^  ^  upon  appealing  to  (i),  (ii)  and  (iii)  of  Lemma 
4.1. 

A  slightly  more  computationally  convenient  form  of  estimator  than  that 

2 

given  by  Theorem  4.5  would  be  obtained  if  the  integrals  involved  in  I  rig 
could  be  replaced  by  discrete  sums. 

If  g  is  a  (sufficiently  smooth)  function  defined  on  an  edge,  y  say, 
of  an  element,  and  if  Qii*»*»Qk  are  the  nodes  which  lie  on  y  then  let  us 
write 

Bg»o  v  *  (lYl  l  lg(Q.)l 2)1/2* 

U,Y  j«l  J 

If  A  is  an  element  and  g  is  defined  on  94  fl  T  then  we  shall  write 

,8,o.aanr  '  (  |  • 

fy  an  edge  of  A') 

1  Ycr  > 

Corollary  4.6.  Theorem  4.5  remains  true  if  E  is  replaced  by  E*  where 

E*  '  (I  l4l(lr,o,aanr)2^2' 

A 

Proof .  Let  r*  be  as  in  (iv)  of  Lemma  4.1.  Since  any  edge  y  can  have  no 
more  than  C  *  C(K)  nodes  on  it  (see  Lemma  1.4.1) 


But 


C_1E*  < 


(I  I A 1 1 r* I 
A 


2  f/2 

o,d&nrJ 


<  CE*. 


|  l4l,t*,o,S4nr  ‘  c  |  l4|f,r,o,«,4nr+ ,r-r*’o,54nr) 

<  C(E2+I  lillr-r*.2i84nr) 

<  C(E2  +  o(h<t>)2) 
by  (iv)  of  Lemma  4.1.  Likewise 

E2  <  C(l  |AHr*82>dAnr+  o(h(P)2)). 

The  corollary  now  follows  by  virtue  of  the  assumed  lower  bound  on  llu-ull^  ^ 
in  the  theorem.  □ 


Since  we  are  only  concerned  here  with  equivalent  estimators  for  llu-ull 


there  are  many  other  modifications  that  can  be  made  to  E  or  E*  without 


affecting  their  equivalent  estimator  property.  Let  us  just  mention  one  other 


such  modification  which  we  shall  say  a  little  more  about  in  Part  II.  Using 


the  same  notation  as  above,  define 


«gi*\  *  (|yl  max  IgCQJl2/2 


f  **  \2  r  **  i2  V?l 

“ax  ('S'o  )  +  “ax  (»g"0  Y)  2 

y  *  '  y  *  ’ 

a  vertical  edge  of  A")  ry  a  horizontal  edge  of  A-j 

ycr  J  ^  yep  > 


i 


igi 


o.aAnr 


V 
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postdoctoral  level,  in  conjunction  with  the  Interdisciplinary  Applied 
Mathematics  Program  and  the  programs  of  the  Mathematics  and  Computer 
Science  Departments.  This  includes  active  collaboration  with  government 
agencies  such  as  the  National  Bureau  of  Standards. 

•  To  be  an  international  center  of  study  and  research  for  foreign 
students  in  numerical  mathematics  who  are  supported  by  foreign 
governments  or  exchange  agencies  (Fulbright,  etc.). 

Further  information  may  be  obtained  from  Professor  I.  BabuSka,  Chairman, 
Laboratory  for  Numerical  Analysis,  Institute  for  Physical  Science  and 
Technology,  University  of  Maryland,  College  Park,  Maryland  20742. 
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